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ABSTRACT 

A new chemical model is presented that simulates fully-coupled gas-phase, grain- 
surface and bulk-ice chemistry in hot cores. Glycine (NH 2 CH 2 COOH), the simplest 
amino acid, and related molecules such as glycinal, propionic acid and propanal, are 
included in the chemical network. Glycine is found to form in moderate abundance 
within and upon dust-grain ices via three radical-addition mechanisms, with no single 
mechanism strongly dominant. Glycine production in the ice occurs over temperatures 
~40-120 K. Peak gas-phase glycine fractional abundances lie in the range 8 x 10 -11 - 
8 x 10" 9 , occuring at ~200 K, the evaporation temperature of glycine. A gas-phase 
mechanism for glycine production is tested and found insignificant, even under op- 
timal conditions. A new spectroscopic radiative-transfer model is used, allowing the 
translation and comparison of the chemical-model results with observations of specific 
sources. Comparison with the nearby hot-core source NGC 6334 IRS 1 shows excellent 
agreement with integrated line intensities of observed species, including methyl for- 
mate. The results for glycine are consistent with the current lack of a detection of this 
molecule toward other sources; the high evaporation temperature of glycine renders 
the emission region extremely compact. Glycine detection with ALMA is predicted 
to be highly plausible, for bright, nearby sources with narrow emission lines. Pho- 
todissociation of water and subsequent hydrogen-abstraction from organic molecules 
by OH, and NH2, are crucial to the build-up of complex organic species in the ice. 
The inclusion of alternative branches within the network of radical-addition reactions 
appears important to the abundances of hot-core molecules; less favorable branching 
ratios may remedy the anomalously high abundance of glycolaldehyde predicted by 
this and previous models. 



Subject headings: Astrochemistry, Astrobiology, ISM: molecules, molecular processes, 
ISM: lines and bands, Radiative transfer 
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1. Introduction 

Millimeter, sub-millimeter, and centimeter wavelength observations of star- forming sources 
known as "hot cores" reveal a range of complex organic molecules, through the identification of 
molecular rotational emission spectra (Herbst & van Dishoeck 2009). Many of the molecules 
detected toward hot cores are highly saturated (i.e. hydrogen-rich), and a significant number also 
contain oxygen and/or nitrogen; examples include methanol (CH 3 OH), formic acid (HCOOH), 
and the ubiquitous methyl formate (HCOOCH 3 ), as well as the recently detected aminoacetonitrile 
(NH 2 CH 2 CN, Belloche et al. 2008), which bears both amine and nitrile groups. The breadth of 
functional groups borne by molecules detected toward hot cores (as well as the large number of 
unidentified emission lines) suggests that other, as-yet undetected species may also be present, 
including biologically-significant molecules such as amino acids - the building blocks of proteins. 

The search for the simplest amino acid, glycine (NH 2 CH 2 COOH), toward interstellar sources 
has not so far yielded a detection (Snyder et al. 2005); recent upper limits toward Sgr B2(N), 
a spectrally- and chemically -rich hot-core source in the Galactic Center, are around 10 14 - 10 15 
cm" 2 for beam-averaged values obtained with the Mopra telescope (Cunningham et al. 2007) and 
the ATCA array (Jones et al. 2007) for glycine conformer I. Upper limits on the higher-energy 
conformer II from the same studies are around 10 13 - 10 14 cm" 2 . However, numerous amino acids 
have been found in meteorites, and glycine was recently detected in a sample taken from comet 
Wild 2 (Elsila et al. 2009); it has indeed been postulated that terrestrial amino acids may have an 
interstellar origin (e.g. Ehrenfreund & Charnley 2000). 

It is likely that amino acids and other biologically-relevant species will, in the coming years, 
provide plausible targets for searches with the ALMA telescope, which can provide much greater 
sensitivity and spatial resolution than was previously available. It is therefore valuable to extend 
existing astrochemical models to simulate the putative formation of glycine and related species 
in hot cores, both to make quantitiative predictions and to assess the viability of observational 
candidate sources. 

Current models of hot-core chemistry place great emphasis on the formation of organic molecules 
within or upon dust-grain ice mantles, as well as in the gas phase. Garrod & Herbst (2006, hereafter 
G&H) showed that typically-observed molecules such as methyl formate, while difficult to form 
by gas-phase mechanisms, could be produced in appropriate abundances through grain-surface 
chemistry at intermediate temperatures, prior to the complete sublimation of the dust-grain ices. 
The radicals HCO and CH 3 were found to become mobile around 20 - 40 K, reacting with other 
radicals associated with formaldehyde and methanol destruction, to produce formic acid, dimethyl 
ether and methyl formate. Cosmic ray-induced photolysis of the molecular ice mantles was found 
to be the dominant source of these functional-group radicals. Garrod, Widicus Weaver & Herbst 
(2008, hereafter GWH) extended the network to include a much larger number of surface radi- 
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cal species, allowing the study of more complex molecules; methanol was found to be the key 
molecule for the supply of the constituent functional groups of many such species. 

Recent experimental work by Oberg et al. (2009a) has provided direct evidence that a range of 
complex organic molecules detected in star-forming regions may be readily produced within UV- 
irradiated organic ices. Other studies focused specifically on amino-acid formation have elucidated 
a number of potential chemical pathways for glycine within laboratory ice mixtures. Experiments 
by Bernstein et al. (2002) and Munoz Caro et al. (2002) showed that UV irradiation of interstellar 
ice analogs containing H 2 0, CH 3 OH, CO, C0 2 , NH 3 and/or HCN could produce significant quan- 
tities of amino acids, including glycine. Sorrell (2001) had suggested a route involving the reaction 
of the NH 2 radical with acetic acid (CH 3 COOH), itself formed through the addition of CH 3 and 
HOCO, while Woon (2002) suggested that the addition reaction between radicals NH 2 CH 2 and 
HOCO could explain the formation of glycine, with these species produced in the ices by succes- 
sive hydrogenation of HCN and by the addition of CO and OH, respectively. Experimental work 
by Holtom et al. (2005) involving the electron bombardment of C0 2 /NH 2 CH 3 ices confirmed the 
viability of Woon's mechanism, while work by Elsila et al. (2007), extending the investigation 
of Bernstein et al. (2002), used isotopic labeling to trace the origins of the constituent atoms in 
the product molecules. Their findings indicated that Woon's mechanism was consistent with the 
origins of the nitrogen and central carbon in the glycine molecule, while the origin of the acidic 
carbon was more consistent with the formation of glycine occurring via a nitrile precursor, i.e. 
amino acetonitrile (NH 2 CH 2 CN). However, a Strecker-type formation route was found to be of 
only minor importance. 

In order to investigate the formation of glycine within dust-grain ices under astrophysical con- 
ditions, the extension of the chemical networks of G&H and GWH - as well as that of Belloche et 
al. (2009), who treated the formation of amino acetonitrile, n-propyl cyanide, and ethyl formate - 
is continued with the addition of glycine, its amino aldehyde glycinal, and a collection of related 
species such as propionic acid and propanal. The new formation mechanisms for these species fol- 
low the general treatments of previous networks, which allowed radicals on the grains, frequently 
produced by UV photodissociation or chemical abstraction processes, to meet via thermal diffusion 
and thus form more complex structures. The mechanism suggested by Woon (2002) is included 
in the new network, as well as several other likely radical-addition mechanisms. Furthermore, a 
gas-phase mechanism for glycine production is tested, by the inclusion of reactions between proto- 
nated hydroxylamine (NH 3 OH + ) and acetic acid (CH 3 COOH), followed by dissociative electronic 
recombination, as suggested by Blagojevic et al. (2003). All new species in the network are af- 
forded both formation and destruction mechanisms on the grains and in the gas-phase. A more 
comprehensive hydrogen-abstraction reaction set is also included in the network. 

To apply this chemical network to hot-core conditions, a new three-phase astrochemical 
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model, MAGICKAL, is used, which includes a fully-active chemistry within the bulk of the ice 
mantle (see Section 2.1). The physico-chemical evolution of the core is traced starting from the 
free-fall collapse of a cloud, using the dust-temperature/visual-extinction relationship given by 
G&P, through the subsequent warm-up of the dense core from 8 - 400 K. 

In order to simulate the emission from glycine and a selection of other molecules, the re- 
sults of these models are mapped to existing observationally-determined temperature and density 
profiles for specific sources. The equation of radiative transfer is then solved numerically, under 
the assumption of local thermodynamic equilibrium, along both on- and off- source lines of sight, 
allowing a determination of line strengths and emission radii. This emission is then convolved 
with telescope beams appropriate to current instruments, allowing the detectability of glycine to 
be assessed. The nearby hot-core source NGC 6334 IRS1 is specifically considered, due both to 
the availability of physical profiles and molecular observations in the literature, and to its spec- 
troscopic suitability. Comparison is also made with the recent non-detection of hydroxylamine 
(NH 2 OH) in another source (Pulliam et al. 2012), with implications for the putative formation of 
glycine in the gas phase. 

Section 2 describes the new chemical and physical models; the results are presented in Section 
3. The results of the inclusion of gas-phase formation mechanisms for glycine are described in 
Section 4. The spectroscopic model and results for various molecules are presented in Section 5. 
A broad discussion of the results is presented in Section 6; a list of general conclusions is provided 
in Section 7. 



2. Chemical modeling 

To model the hot-core chemistry, a new astrochemical code is employed, named MAGICKAL 
(Model for Astrophysical Gas and Ice Chemical Kinetics And Layering). This code uses a rate- 
equation/modified rate-equation approach to solve the coupled gas-phase, dust grain-surface and 
ice-mantle chemistry. The essential functionality of previous approaches (Garrod & Pauly, 2011) 
is retained, while a number of additional features are added, as outlined below. 



2.1. Three-phase model 

The three-phase model used here is described in detail by Garrod & Pauly (201 1), and is based 
on the approach of Hasegawa & Herbst (1993). In contrast to both of those methods, MAGICKAL 
simulates an active ice-mantle chemistry in addition to that occurring in the surface layer of the 
ice or on the dust-grain surface. Previous two-phase hot-core models, such as those of G&H, 
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GWH, and Belloche et al. (2009), have implicitly assumed a combined surface/mantle chemistry, 
using surface kinetics to simulate all dust-grain chemical processes, including those within the bulk 
ice. Although satisfactory results may be obtained through this simplification, it has become clear 
through experimental means (e.g. Oberg et al., 2009a) that the promotion of complex-molecule 
formation within organic ices by UV photolysis involves sub-surface processes, even in relatively 
thin ices (on the order of 10 monolayers). 

In order to simulate chemical kinetics within the ice mantles, it is here assumed that the mo- 
bility of chemical reactants in the bulk ice occurs through a barrier-mediated thermal diffusion 
process. Other authors have assumed a swapping mechanism to be active in laboratory ice experi- 
ments (Oberg et al. 2009; Fayolle et al. 201 1), which is also implicitly assumed in this model. The 
swapping rate is parameterized in the same way as the surface thermal hopping, i.e. 

*swap(0 = Mi) exp[-£ swap (/)/r rf ] ( 1 ) 

where u is the characteristic vibrational frequency of an harmonic oscillator, T d is the dust tem- 
perature, and inswap is an energy barrier associated with the swapping of species i with an adjacent 
water molecule, in line with the adoption of surface binding energies appropriate to amorphous 
water ice, following G&H and GWH. 

Ratios between the surface diffusion barrier and the desorption energy of any particular 
species bound to an ice or any other surface are poorly constrained, and various authors have 
used values in the approximate range of 0.3 - 0.8. Garrod & Pauly (201 1) found that the produc- 
tion of interstellar CO2 ice could be well reproduced by the adoption of ratios E^-.E^ < 0.4; here 
a value of 0.35 is used. Thence, a value of E swap :E des = 0.7 is assumed for the swapping barrier, on 
the simplistic consideration that a bulk-ice molecule is bound to approximately twice the number 
of binding partners as a surface molecule. In analogy to the surface chemical reaction rates, the 
rate of reaction between mantle-species A and B is given by: 

ab =N m (A) N m (B) [fc swap (A) + fc swap (fl)] /N M (2) 

where N m (i) is the population of species i in the ice mantle, and N M is the total population of the 
ice mantle, including all species. 

As discussed by Garrod & Pauly (2011), in the three-phase model, material in the surface 
layer is transferred to the bulk ice according to the net rate of deposition of material onto the grain 
from the gas phase, with a similar, alternative transfer occurring in reverse in the case of a net 
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desorption of mantle species into the gas phase. These transfers represent the covering or exposure 
of ice material as the boundary of the ice-surface layer is constantly re-defined within the model, 
rather than the physical transport of material. However, the inclusion of bulk diffusion within the 
ice should allow a true exchange of chemical species between surface and mantle, through the 
aforementioned swapping mechanism. 

This surface-mantle swapping is a pair-wise process, and therefore requires the construction 
of rates for each chemical species in the mantle and in the surface layer, which must collectively 
produce no net transfer between the two phases. This zero net rate would be assured using an 
explicit pair-wise coupling of every species in the surface with every species in the grain mantle; 
however, such an approach is found to be computationally exhaustive, and requires estimates of 
the pair-wise positions of the swapped species that a rate-based model cannot provide, other than 
through purely statistical considerations. 

Instead, the swapping rates from mantle to surface are calculated first, with rates analogous to 
that of equation (2). The total rate in this direction is then matched with an identical total transfer 
rate from surface to mantle, apportioned to all the surface species according to their fractional 
contribution to the total surface-layer population. Thus: 

N s 

^swap.mO) ~ N m (f) — &swap(0 (3) 



^swap,s(0 — jy ' / v ^swap,mQ) (4) 
^ all j 

where i? SW ap,mO) is the rate of mantle-to-surface swapping of mantle species /, 7? swap s (z) is the 
rate of surface-to-mantle swapping of surface species i, N m {i) and N s (i) are the mantle and surface 
populations, respectively, of species /, and N M and N s are the total mantle and surface populations. 
The ratio N S /N M in equation (3) is taken as unity when it exceeds this value, representing the case 
where some surface species reside on the dust-grain surface itself and have no underlying mantle. 

Because the model does not assign specific positions to each atom or molecule in the mantle, 
the rate of mantle-to-surface swapping, ^ swap , m (z), is required to be representative not simply of the 
final diffusion step, but of the entire chain of diffusion steps by which a mantle species at arbitrary 
distance from the surface layer will reach the surface. The use of /c swap in eq. (3) as the rate of a 
typical, or average, step in this process is accurate so long as the final diffusion step from mantle 
to surface is relatively fast, i.e. that intra-mantle diffusion is the rate-limiting step in the entire 
journey from mantle to surface. Under the assumption that the ice is predominantly water ice (see 
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above), the final surface-to-mantle step will be fast for species whose binding energies are less 
than or similar to that of water. For the small minority of species with significantly greater binding 
energies than this, for which there would thus be an associated endothermicity to the mantle- 
surface swap with a water molecule (see e.g. Oberg et al. 2009b), the large absolute size of the 
swapping barrier would render the rate of mantle-to-surface diffusion negligible compared to other 
processes, such as the mantle-to- surface transfer associated with net surface desorption. 

These rates defined in eqs. (3) and (4) are also considered in the loss and gain rates used for 
the modified-rate calculations of Garrod (2008), which are applied in this model (using "method 
C"). Garrod et al. (2009) showed that this method produces a very close approximation to the 
results of exact Monte Carlo methods, using a large chemical network under a range of physical 
conditions. 



2.2. Chemistry 

The chemical network used in this model is based on that of GWH and the subsequent addi- 
tions of Belloche et al. (2009), which featured a treatment for nitriles, including amino acetonitrile 
(NH 2 CH 2 CN). These networks assumed grain-surface radical-radical formation mechanisms for 
newly-added complex organic molecules, as well as destruction mechanisms such as grain-surface 
and gas-phase photodissociation (by both the ambient and cosmic ray-induced UV fields), gas- 
phase ion-molecule reactions (including protonation) and the related electronic recombinations of 
the products, as well as accretion from the gas-phase onto the grain surfaces, and thermal and 
n on-thermal desorpt ion from the grains. The initial elemental abundances used are shown in Table 



[B lLaas et al.l (|201l|) constructed a network of reactions to treat the formation of frans-methyl for- 
mate, the higher energy conformer of the commonly observed czs-methyl formate; however, that 
study found no significant effect on the abundances of species other than trans-MF, so the omission 
of the related chemistry is not expected to have any effect on the models presented here. 

Following the same methodology, formation and destruction processes are included for glycine 
and a suite of related molecules, including glycinal, propanal and propionic acid. Radical-radical 
reactions for a selection of these species are shown in Table [2] Whereas in the previous models, 
a single complex molecule was typically assumed to be the only possible product of the addition 
of key grain-surface radicals, here alternative product branches are included, for the case where 
such a reaction could also result in a pair of stable products, and where the reaction is also exother- 
mic. Thus, for the reaction HCO + CH 3 ->■ HCOOCH 3 , which was found by GWH to be the 
predominant formation route for methyl formate, alternative product branches of [CO + CH 3 OH] 
and [H 2 CO + H 2 CO] are now also included. The branching ratios are in general poorly defined, 
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and are likely to be significantly different between the gas- and solid-phase processes; statistical 
ratios are therefore adopted (e.g. 1:1:1 for HCO + CH 3 0). In addition, the same such two-product 
reactions are now also included in the gas-phase chemistry, using measured ratesQ where available; 
otherwise, a uniform second-order rate coefficient of 10" 11 cm 3 s" 1 is assumed. Such processes 
may become important in the late, high-temperature stage of hot-core evolution, during which the 
evaporation, protonation and subsequent electronic recombination may lead to significant quanti- 
ties of molecular radicals being present in the gas phase. However, the single-product branches 
are omitted from the gas-phase chemistry, because of their negligible rates, due to the absence of a 
grain surface to carry away excess energy and thereby stabilize the product. 

Following previous models, grain-surface hydrogen-abstraction reactions are added for reac- 
tions between the new complex species and the most significant grain-surface radicals such as H, 
OH, CH 3 , CH 2 OH, CH 3 and NH 2 . Table [3] shows activation-energy barrier information used in 
the model, for a selection of (primarily hydrogen- abstract ion) reactions occurring on the grains. 
For atomic -hydrogen reactions relating to methanol formation, for which more complex fits to 
theoretical rates are employed, the simpler rectangular-barrier parameter fits are quoted, as for all 
other barrier-mediated reactions. The method of Hasegawa et al. (1992) is used to calculate the 
tunneling and thermal rates associated with these parameters, adopting the faster of the two for any 
given temperature. For the majority of reactions, only gas-phase data are available, in which case 
the activation energy derived from fits to the lowest-temperature data-sets available are adopted 
for the surface/mantle reactions, assuming a barrier width of 1 A. In cases where no data exist, 
the activation energy, E A , is estimated using the Evans-Polanyi relationship, which states that the 
difference in E A between reactions of a similar type is proportional to the difference between the 
enthalpy changes associated with each reaction. Parameters for the Evans-Polanyi relationship are 
obtained from Dean & Bozzelli (1999), or from fits using similar reactions, indicated in Table 
[3] Where the enthalpy change of reaction is unknown, the activation energy of a similar reac- 
tion is used directly. There is evidence that the presence of a surrounding water-ice matrix may 
alter activation energy barriers for specific product branches (e.g. Goumans 2011, in the case of 
formaldehyde-related chemistry), or lower the internal barriers to the structural re-arrangement of 
a product (e.g. Duvernay et al. 2005, in the case of the isomerization of cyanamide). It is there- 
fore plausible that some of the gas-phase activation energies used in this model over-estimate the 
barriers in the solid phase; however, such effects are difficult to quantify without a specific com- 
putational or experimental study of each solid-phase reaction and its low-temperature tunneling 
characteristics. 

For self-consistency, all reactions of surface radicals with complex species are now assumed 
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to result in hydrogen abstraction, in contrast to GWH, who included a small selection of alternative 
substitution reactions for aldehyde group-bearing species. These reactions typically provided only 
minor destruction or formation routes for the species involved, and their removal has no significant 
influence over the glycine chemistry described below. 

All processes related to grain-surface chemistry outlined above are likewise incorporated into 
the chemistry of the bulk ice, with identical activation energies (where present), while their reaction 
kinetics are instead defined by the bulk-ice diffusion parameters discussed in Section 2.1. Photodis- 
sociation of ice mantles is also included to all depths, without attenuation, in line with previous 
models. This is a reasonable first-order approximation; using the photon absorption cross-section 
of amorphous water ice for the first absorption band measured by Mason et al. (2006), Andersson 
& van Dishoeck (2008) calculate an absorption probability per photon of 0.007 per monolayer of 
water ice. The maximum quantity of water ice formed on the grains in this model corresponds to 
~1 14 monolayers; this would result in a UV field strength at the deepest ice layer of ~45% of the 
value at the surface. This corresponds to a mean value of 69% of the assumed field strength, for 
the ice mantle as a whole. However, the ices formed on interstellar dust-grain surfaces constitute 
a mixture of various molecular species, and the absorption cross-sections for many of these are 
currently unknown or poorly defined. 

No desorption processes are allowed from the bulk ice; thus desorption of these species must 
proceed via a prior transfer of material from mantle to ice surface. Reaction-diffusion competition 
is also included for reactions mediated by activation-energy barriers (following G&P), for both 
surface and bulk-ice reactions; this is perhaps the largest departure from the methods of GWH, as 
this change can produce reaction efficiencies many orders of magnitude different from the simple 
Arrhenius treatment of previous models. 

Surface desorption processes included are thermal evaporation (as defined by the binding/desorption 
energy, E des ), reactive desorption as described by Garrod et al. (2007), and photo-desorption us- 
ing measured rates determined by Oberg et al. (2009c,d); following G&P, all surface species are 
assumed to photodesorb with an efficiency of 10" 3 , in the case where values are not provided by 
those authors. 

The indirect mechanism for C0 2 formation proposed by G&P is also included for reactions 
both on the surfaces and within the ice mantles. 

The reproduction of the high observed abundance (~10~ 8 with respect to H 2 ) of methyl for- 
mate (HCOOCH3) in hot cores has proved somewhat problematic in previous astrochemical mod- 
els. Garrod & Herbst (2006) showed that grain-surface formation mechanisms could produce 
adequate quantities of methyl formate (MF) to explain observations; however, in some cases it 
was found (Garrod et al. 2008; Belloche et al. 2009) that the evaporation of MF at relatively low 
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temperatures could result (dependent on the warm-up timescale) in only a small fraction of this 
material surviving the transition to the gas-phase. Belloche et al. (2009) adopted a higher binding 
energy for MF (5200 K vs. 4100 K), citing the importance of hydrogen bonding in the estimation 
of this value; they found that this value produced results appropriate to observational abundance 
values for methyl formate. This value is also adopted here. Table @] shows a selection of other 
binding energies used in the present model. 



2.3. Methanol chemistry 

As found by GWH and Laas et al. (2011), the rates and branching ratios of methanol pho- 
todissociation in the ice mantles are crucial to the production rates and ratios of various complex 
organic molecules detected in hot cores, as the ice-mantle methanol provides a large supply of 
structurally-complex organic material. 

The present model uses the methanol branching ratios derived from the experiments of Oberg 
et al. (2009a) for products [CH 2 OH+H]:[CH 3 0+H]:[CH3+OH] of 5:1:1. In the absence of other 
information, these values are applied to both the direct UV photodissociation by the interstellar 
radiation field, and to photodissociation by the secondary UV photons produced by cosmic rays, 
the latter of which is most important at the high visual extinctions of the hot-core phase. These 
ratios are applied throughout to the ice and gas phases, although, as found by Laas et al. (201 1), the 
importance of gas-phase methanol photodissociation to complex molecule production is minimal. 

Activation-energy treatments for CO and H 2 CO use best-fit expressions to uni-molecular 
and bi-molecular rates recently acquired through quantum transition state theory (Goumans, priv. 
comm.; Andersson, priv. comm.). Hydrogen abstraction from H 2 CO and CH 3 OH are included, 
mediated by activation energy barriers, while hydrogen abstraction from HCO and CH 3 0/CH 2 OH 
are also allowed (producing H 2 + CO/H 2 CO), assuming branching ratios of 1:1 with the associated 
hydrogen-addition channels (i.e. no activation energy). 



2.4. Glycine chemistry 

2.4.1. Grain-surface/mantle chemistry 

A number of potential routes have been suggested in the literature for the formation of glycine 
in interstellar ices, based on laboratory evidence or theoretical considerations. Experiments by 
Bernstein et al. (2002) and Munoz Caro et al. (2002) showed that UV irradiation of interstellar ice 
analogs containing H 2 0, CH 3 OH, CO, C0 2 , NH 3 and/or HCN could produce significant quantities 
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of amino acids, including glycine. Sorrell (2001) had suggested a route involving the reaction of the 
NH 2 radical with acetic acid (CH 3 COOH), itself formed through the addition of CH 3 and HOCO, 
while Woon (2002) suggested that the addition reaction between radicals NH 2 CH 2 and HOCO 
could explain the formation of glycine, with these species produced in the ices by successive 
hydrogenation of HCN and by the addition of CO to OH, respectively. Experimental work by 
Holtom et al. (2005) involving the electron bombardment of C0 2 /NH 2 CH 3 ices confirmed the 
viability of Woon's mechanism. Elsila et al. (2007) found Woon's mechanism to be consistent 
with the origins of the nitrogen and central carbon in the glycine molecule, while the origin of the 
acidic carbon was more consistent with the formation of glycine occurring via a nitrile precursor, 
i.e. amino acetonitrile (NH 2 CH 2 CN). 

Here, we include a number of new radical-addition reactions relating to glycine and similar 
species. Each hydrogenation step suggested by Woon (2002) is included, although all but the final 
addition of NH 2 CH 2 + HOCO was already present in the network. However, other routes are 
also included, such as the addition reaction NH 2 CH 2 + HCO — > NH 2 CH 2 CHO (glycinal), whose 
product may (for example) be photodissociated to produce the radical NH 2 CH 2 CO, which may 
further react with OH to produce glycine. There are four final processes that lead to glycine in this 
network, 



although there are several ways in which each radical may come to be formed, either through 
chemical reaction or by photodissociation of a stable molecule. Reaction (5) is also complemented 
by exothermic radical-radical reactions with species (e.g. HCO) that may donate a hydrogen atom 
to form two stable molecules, as described in section 2.2. Reaction (5) typically follows the pho- 
todissociation of, or H-abstraction from, glycine itself; reactions (6)-(8) are thus the major potential 
formation routes. 

Woon suggests the formation of HOCO by the reaction of CO and OH. This process is difficult 
to model accurately, due to the complex nature of the transition states; the reaction may also lead 
to the formation of C0 2 + H (see G&R and references therein), or may dissociate back to CO + 
OH, dependent on the rapidity of energy loss to the surface. As such, the OH + CO reaction is 
modeled separately for each set of products, assuming an activation energy barrier of 2500 K for 
the HOCO branch, following the approach of GWH. In practice, the formation of HOCO in the 
model is dominated by the photodissociation of or chemical hydrogen abstraction from formic acid 
(HCOOH) in the ice, which is formed predominantly by HCO + OH additions (see GWH). 



H + NHCH 2 COOH -> NH 2 CH 2 COOH 
NH 2 + CH 2 COOH NH 2 CH 2 COOH 
NH 2 CH 2 + HOCO ->• NH 2 CH 2 COOH 
NH 2 CH 2 CO + OH ->• NH 2 CH 2 COOH 



(5) 
(6) 
(7) 
(8) 
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Also included in the model is a selection of species similar to glycine, but with alternative 
functional groups. For example, glycinal (NH 2 CH 2 CHO) and propionic acid (C 2 H 5 COOH) are 
modeled, through analogous reactions to those of glycine, while new dust-grain formation routes 
are added for propionaldehyde (or propanal, C 2 H 5 CHO), which has been detected toward Sgr 
B2(N) (Hollis et al. 2004). The inclusion of these other species ensures that the importance of 
reactions leading to glycine formation is not over-estimated due to the omission of competing 
reaction pathways. 



2.4.2. Gas-phase chemistry 

As for other complex species, the chemical network includes gas-phase destruction routes 
for glycine that consist of UV photodissociation, either by the cosmic ray-induced or standard 
interstellar radiation field, or by reactions with atomic or simple molecular ions. Table |5] shows the 
newly-included reactions for glycine and its related gas-phase products. Electronic recombination 
of protonated glycine results in the re-formation of glycine with an efficiency of 5 %, but all the 
glycine in the gas-phase is originally formed on the grains. Consideration of a potential gas-phase 
formation route is given in the models of Section 5. 

2.5. Physical model 

The physical model used here broadly follows that of G&H and later papers, in which a 
collapse phase is followed by a static warm-up phase. In the first phase, the density increases from 
«h = 3x 10 3 to 10 7 cm" 3 , under the free-fall collapse expression given by Nejad et al. (1990). 
However, in contrast to previous hot-core models, while the gas temperature is held at a constant 
10 K, the temperature of the dust varies with visual extinction, according to equation (17) of G&P, 
which describes the heating of the dust by the interstellar radiation field, assuming a dust-grain 
radius of 0.1 fj,m. At the initial visual extinction of A v = 2, the dust temperature is ~16 K; during 
collapse, this falls to a minimum temperature of 8 K. 

During the subsequent warm-up phase, the core is heated from 8 - 400 K, with density held 
at 10 7 cm" 3 . Three warm-up timescales are adopted, which are the same as previous models, in 
so far as they reach 200 K at the same point in time (Table©. However, the warm-up is extended 
beyond this time, until 400 K is reached, so that the later evaporation and gas-phase behavior of 
more strongly-bound species, such as glycine, may be studied over a longer period. During the 
warm-up phase, gas and dust temperatures are assumed to be well coupled due to the high density; 
the gas kinetic temperature follows that of the dust grains when r dust >10 K. Table [6] shows the 
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details of the three warm-up phase physical models. 

Table [7] shows the composition of the dust-grain ices at the end of the collapse phase (corre- 
sponding to the starting values for each of the following warm-up models), as well as observational 
values for comparison. 

3. Results 

Figures 1 , 2 and 3 show the behavior of fractional abundances with respect to total hydrogen 
(«h) of a selection of simple and complex molecules, for the fast, medium and slow warm-up mod- 
els, respectively. Peak gas-phase abundances are shown in Table [8} along with the corresponding 
system temperatures for each peak value, which frequently correspond to the temperature of peak 
evaporation for species formed on the dust grains. 

In Figs. 1-3, panel (a) shows some of the main ice constituents and a representative ion, 
HCO + . Panels (b) and (c) show some complex organics typically observed toward hot cores. Panel 
(d) shows some commonly observed nitrogen-bearing species as well as the as-yet undetected 
hydroxylamine (NH 2 OH). Panel (e) shows a selection of complex species containing a carbonyl 
group (C=0), including acetone, (CH 3 ) 2 CO. Panel (f) shows complex species derived from the 
larger photo-fragments of methanol. Panel (g) shows a selection of complex organics bearing 
an amine group (NH 2 ), while panel (h) shows glycine itself (NH 2 CH 2 COOH) along with several 
similar or related molecules. Not all species shown in Table [8]are plotted in Figs. 1-3. 

The results for a selection of important species and their progenitors are discussed in detail 
below. Due to the large number of molecules included in the model, the results for many molecules 
are presented in Figs. 1-3 and Table [8] without comment. For a description of the most important 
reaction mechanisms relating to these species, the reader is referred to the analysis presented by 
Garrod et al. (2008); however, some general differences between the current and previous models 
are discussed in section 6. 



3.1. Glycine and related species 

3.1.1. Glycinal 

Glycinal (NH 2 CH 2 CHO) is the amino aldehyde corresponding to the amino acid glycine. Its 
formation in these chemical models, shown in panel (h) of Figs. 1-3, begins at around 25 K, when 
the radical HCO becomes mobile and reacts with CH 2 NH 2 . The latter species is formed by the 
hydrogenation of CH 2 NH, which builds up at around 15 K via the addition of NH and CH 2 radicals. 
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The HCO radical is mainly formed by the abstraction of hydrogen from formaldehyde (H 2 CO) by 
OH, which is a barrierless process; this occurs until formaldehyde evaporates at around 40 K (see 
G&H, GWH, and section 3.2, below). A later and somewhat stronger glycinal formation route 
operates at ~40 K, when NH 2 radicals add to CH 2 CHO radicals formed by the photodissociation 
of acetaldehyde (CH 3 CHO). 

Glycinal formation thus occurs at rather lower temperatures than that of glycine (see below). 
However, both are at least partially dependent on the earlier formation of acetaldehyde at lower 
temperatures to provide the CH 2 CO structure present in both. Peak gas-phase glycinal abundances 
are a few 10" 9 to 10" 8 with respect to total hydrogen; it is thus of similar abundance to the more 
prevalent observed organic molecules in hot cores. No mm/sub-mm spectral information is appar- 
ently available for glycinal at present. 



3.1.2. Glycine 

Glycine formed on the grains is seen to be released into the gas phase at a temperature of 
around 200 K, attaining a peak gas-phase abundance of 8.6 x 10 -11 , 4.6 x 10" 10 , and 8.1 x 10" 9 
with respect to total hydrogen (see Tableland Figs. 1-3 panel h), for the fast, medium, and slow 
warm-up timescales, respectively. 

The three main production routes included for glycine (Eqs. 6 - 8) are all active, and all 
contribute a proportion of the glycine later released into the gas phase. Significant dust-grain 
surface/mantle formation of glycine begins at temperatures of 40 - 50 K, and may continue up to 
around 120 K. 

For the medium warm-up timescale case, at around 40 K the OH radicals in the ice, produced 
by the photodissociation of water molecules, are sufficiently mobile to find glycinal (NH 2 CH 2 CHO), 
which is abundant at around 10" 9 - 10" 8 with respect to total hydrogen. The OH radical easily ab- 
stracts a hydrogen atom from glycinal; the activation energy barrier is sufficiently low (591 K) 
compared to the diffusion barrier for OH that the reaction occurs before OH has an opportunity to 
diffuse away from its reaction partner. (This effect indeed confers an efficiency close to unity on 
most such OH-initiated hydrogen-abstraction reactions). The resulting radical, NH 2 CH 2 CO, then 
reacts with another OH to produce glycine (Eq. 8). 

At temperatures of 55 K and above, much greater quantities of glycine are formed; hydrogen 
abstraction from acetic acid (CH 3 COOH) by NH 2 radicals within the ice mantles and the sub- 
sequent addition of the resultant CH 2 COOH radicals to NH 2 produces glycine (Eq. 6). The NH 2 
radicals are primarily formed by hydrogen abstraction from ammonia (NH 3 ) by OH radicals, which 
are themselves formed by the photodissociation of water ice. The increased importance to glycine 
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formation of the acetic acid-related path at these temperatures is caused by the greater abundance 
of NH 2 radicals that results from the evaporation of HNCO molecules from the ices, which up until 
that time represent a major reaction partner for NH 2 . The release of HNCO into the gas phase may 
be seen in Figs. 1-3, panel (d). However, over the 55 - 75 K range, Eq. (8) still dominates Eq. (6) 
by a factor of 2 - 3, as the greater availability of NH 2 allows the NH 2 CH 2 CO radical to be formed 
by hydrogen abstraction from glycinal by NH 2 as well as OH. 

At around 75 K, an ice surface (as opposed to mantle) process becomes the strongest con- 
tributor to glycine formation; hydrogen abstraction from CH 3 NH 2 (methylamine) by OH radicals 
produces CH 2 NH 2 , while abstraction from HCOOH (formic acid) by OH produces HOCO. These 
react to produce glycine (Eq. 7) on the ice surface, by the route suggested by Woon (2002). This 
mechanism then competes roughly equally with the reaction of CH 2 COOH and NH 2 on the ice 
surface (Eq. 6) to be the main producer of glycine, up to ~90 K. By this point, most glycine 
formation is taking place on the surfaces of the ices, rather than within the ice mantles. This is 
caused by the increasing rate of mantle-to- surface diffusion among less strongly-bound species, 
including NH 2 , as the temperature rises, coupled with their commensurately faster reaction rates 
on the ice surface. The more strongly-bound products of these reactions, such as glycine itself, 
are "swapped" back into the deeper mantle following their formation, as they are exchanged with 
more mobile, volatile species. 

Between around 90 - 120 K, the reaction of NH 2 and CH 2 COOH (Eq. 6) again takes over as 
the dominant route to glycine on the grains, as CH 3 NH 2 is gradually lost to the gas phase, rendering 
Eq. (7) far less effective. 

The majority of the glycine ultimately released into the gas phase is formed during the 90 - 
100 K period. During this time, Eq. (6) dominates Eq. (7) by a factor of ~2. 

At around 120 K (in this and all models), ammonia (NH 3 ), formic acid (HCOOH), and acetic 
acid (CH3COOH) evaporate fully from the grains, closing down the main glycine formation mech- 
anisms. From this point on, glycine begins to migrate from the bulk ice mantle into the surface, 
while suffering a moderate degree of attrition by hydrogen abstraction and photodissociation - the 
products of which typically evaporate. 

For the fast warm-up model, the temperature dependence of the formation mechanisms is 
broadly similar; however, the importance of each route is different, due to the differing quantities 
of precursor molecules available. Most importantly, the abundance of acetic acid (CH 3 COOH), 
shown in panel (e), is around one order of magnitude lower than in the medium timescale model. As 
found by Garrod et al. (2008), in these models acetic acid is formed predominantly by the addition 
of OH to CH3CO, the product of the abstraction of hydrogen from acetaldehyde (CH 3 CHO) by OH 
itself. The fast warm-up model allows a shorter time period for this process, before acetaldehyde 
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evaporates from the grains, resulting in the lower acetic acid abundance. This in turn renders the 
production of glycine via Eq. (6) generally a factor of a few slower than Eqs. (7) & (8), which, 
overall, provide approximately equal amounts of glycine production in this model. 

For the slow warm-up model, the behavior of glycine production is much more similar to that 
of the medium case; however, here the NH 2 + CH 2 COOH route is strongly dominant over the 60 
- 120 K range during which glycine production is important. This is caused by the greater time 
period for the production of acetic acid. Production via Eq. (8) is also less significant than in the 
two shorter- timescale models, due to the lower abundance of glycinal in the ices. 

In all models, the gas-phase destruction of glycine is entirely dominated by protonation and 
dissociative electronic recombination, as is the case for the majority of complex organics (see 
also GWH); photodissociation has only a small effect, as the only UV radiation field capable of 
influencing the chemistry during the warm-up phase is that induced by cosmic rays. 

These models produce a low to moderate abundance of glycine in the gas phase, in comparison 
with other complex species produced in the models. The peak gas-phase abundance of glycine 
is also very similar to that of acetic acid in each model, albeit with a different temperature of 
evaporation. The question of whether this glycine may be detectable is addressed in Section 5. 



3.1.3. Amino acetonitrile 

Amino acetonitrile (NH 2 CH 2 CN), or AAN, shown in panel (h), has been suggested as a pre- 
cursor to glycine (Elsila et al. 2007). It is formed in these models with a peak gas-phase fractional 
abundance of a few 10" 9 to 10" 8 , consistently for each timescale model. This is around a factor of 
two higher than the values produced in the model of Belloche et al. (2009), and around a factor of 
10 greater than the value derived from observations toward Sgr B2(N) by Belloche et al. (2008). 

In these models, AAN is formed primarily by the addition of NH 2 to the radical CH 2 CN, 
which itself derives from hydrogen abstraction from CH 3 CN molecules by NH 2 at around 60 - 80 
K, producing an abundance of around 10" 9 on the grains. Later, at around 90 K, the evaporation 
of CH 3 CN from the grains results in its gas-phase protonation and recombination with electrons, 
producing CH 2 CN with an efficiency of 50 %. This radical re-accretes onto the grains, to react 
with NH 2 , which produces an AAN fractional abundance of around 10" 8 on the grains. However, 
in light of recent measurements of methanol recombination branching ratios (e.g. Geppert et al. 
2006), it is unclear whether such a high efficiency of CH 2 CN production by gas-phase electronic 
recombination is realistic. A value of 5 %, as adopted in these models for the recombination 
of more complex species, would render the earlier formation of AAN by purely grain- surface 
chemistry the dominant mechanism, and would produce abundances close to the values observed 
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toward Sgr B2(N). 



3.1.4. Propionaldehyde and propionic acid 

Propionaldehyde (or propanal, C 2 H 5 CHO) and propionic acid (or propanoic acid, C 2 H 5 COOH) 
are structurally similar to glycinal and glycine, and neither have been treated in chemical models 
previously. Propionaldehyde is formed in great abundance at low temperatures, by the addition of 
HCO and C2H5 radicals on the dust-grain ice surfaces; abundances as high as 10" 7 with respect to 
total hydrogen are formed in the slow warm-up model. Production is most rapid at around 30 K, 
when the evaporation of grain-surface methane (CH 4 , panel a) is fastest. This results in a signifi- 
cant increase in the abundance of gas-phase acetylene (C 2 H 2 ), which accretes onto the dust grains, 
where it is rapidly hydrogenated to C 2 H 5 , thence reacting to form propionaldehyde. Formation on 
the ice surface at this time is sufficiently rapid that so-called reactive desorption, whereby the en- 
ergy released from the exothermic formation of a surface molecule results in desorption (e.g. Gar- 
rod et al. 2007), produces a gas-phase peak in the abundance of propionaldehyde of approximately 
10" 12 n u . This period of carbon-rich gas-phase chemistry associated with methane evaporation is 
often labeled warm carbon-chain chemistry (WCCC), and has been identified observationally by 
Sakai et al. (2008); however, it is unlikely that the enhancement in propionaldehyde seen in the 
models would be observable at the low abundances predicted by these models. The large quantity 
that remains on the grain surfaces later evaporates at around 120 - 130 K, producing gas-phase 
peak fractional abundances of around 10" 8 . This molecule was detected by Hollis et al. (2004) 
toward Sgr B2(N); however, no specific abundance was determined. 

Propionic acid (C 2 H 5 COOH) is formed on the grains much later than propionaldehyde; most 
production occurs from around 70 K onwards, by the addition of OH and C 2 H 5 CO radicals on the 
ice surface. The latter species is formed primarily through hydrogen abstraction from propionalde- 
hyde by NH 2 in the ice. Thus, propionic acid production appears to be directly dependent on the 
prior formation of its associated aldehyde. 



3.2. Other complex molecules 

3.2.1. Methanol and formaldehyde 

Methanol (CH 3 OH) and formaldehyde (H 2 CO), shown in Figs. 1-3 panel (b), are both present 
in the dust-grain ice mantles in large quantities, following their formation at low temperature during 
the collapse phase (see Table[7]). Formaldehyde remains the primary source of HCO radicals in the 
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ices, due to hydrogen abstraction by OH, until it evaporates at around 40 K, producing a significant 
gas-phase peak - see G&H and GWH. This feeds the gas-phase formation of other species, most 
notably formic acid and methanol. 

Methanol abundance in the ice is essentially static until it eventually evaporates at around 120 
K. But methanol nevertheless provides the majority of organic structure for the formation of more 
complex molecules, in the form of methyl (CH 3 ), methoxy (CH 3 0) and hydroxymethyl groups 
(CH 2 OH). These are formed primarily by photodissociation caused by the cosmic ray-induced 
UV field that penetrates the hot core. Above ~45 K, OH radicals become mobile, and hydrogen 
abstraction to form CH 2 OH becomes the dominant destruction mechanism for methanol, producing 
around twice as much CH 2 OH as the photodissociation route. 

The abundance of methanol in the ices is little affected by the warm-up timescales of each 
model. The fall in the gas-phase abundances is more significant, but the apparent differences be- 
tween each timescale model are mainly due to the longer timescales over which the chemistry 
of the slower warm-up models is simulated, rather than to a difference in the rate of destruc- 
tion. The temperature-abundance relationships may, however, be important to the interpretation of 
observationally-determined abundances and rotational temperatures. The apparently precipitous 
fall in the gas-phase abundances of various other molecules in the slow warm-up model is also 
largely attributable to the longer timescales adopted in that model. 

The eventual evaporation of methanol also indirectly produces a high-temperature peak in 
formaldehyde abundance; the reaction of methanol with protonated methanol produces protonated 
dimethyl ether (CH 3 OCH 4 + ), one of whose electronic recombination products is assumed to be 
formaldehyde in this model. The branching products of this process are not well constrained, pace 
Hamberg et al. (2010), who indeed do not identify H 2 CO as a likely product. The secondary peak 
in formaldehyde abundance may at least be considered unreliable, due to these uncertainties. 



3.2.2. Methyl formate and glycolaldehyde 

The peak abundance of methyl formate (MF, HCOOCH 3 ), shown in Figs. 1-3 panel (b), is 
in line with that expected from observations (around 10" 8 n H ), except for the slow warm-up model. 
The evaporation of methyl formate prior to that of methanol and various other species leaves it 
vulnerable to destructive reactions with molecular ions, whose abundance would otherwise be 
lessened by the presence of other reaction partners. The long timescale between methyl formate 
and methanol evaporation in the slow model therefore leads to significant methyl formate destruc- 
tion; a similar result was noted by GWH and Belloche et al. (2009), in the case of a lower assumed 
MF binding energy. 
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Methyl formate is formed entirely on the dust grains in this model, by the reaction HCO + 
CH 3 0. Production is mainly within the ice mantles, peaking at 25 K, although significant produc- 
tion also occurs on the ice surface. 

Glycolaldehyde (CH 2 (OH)CHO), shown in Figs. 1-3 panel (c), is formed similarly to its a 
structural isomer methyl formate, via HCO + CH 2 OH. Its gas-phase and ice abundance is greater 
than that of MF, due to the more rapid production of CH 2 OH over CH 3 - by a factor of 5 (for 
photodissociation) in this model. Observations of protostellar envelopes, however, do not demon- 
strate this large abundance of glycolaldehyde; the recent detection by J0rgensen et al. (2012) 
toward the Class protostellar binary IRAS 16293-2422, using the ALMA telescope, indicated a 
HCOOCH 3 :CH 2 (OH)CHO ratio of 10 - 15 (although the observational glycolaldehyde rotational 
temperature of 200 - 300 K is consistent with the present model results). 

The discrepancy may be related to the assumed branching ratios of the reactions between 
HCO and either CH 3 or CH 2 OH: 



and 



HCO + CH3O 



HCO + CH 2 OH 



->• HCOOCH3 
^CO + CH 3 OH 
^H 2 CO + H 2 CO 



-»■ CH 2 (OH)CHO 
^CO + CH 3 OH 
^H 2 CO + H 2 CO 



(-423 kJ/mol) 
(-372 kJ/mol) 
(-292kJ/mol) 



(-290 kJ/mol) 
(-276 kJ/mol) 
(-197kJ/mol) 



(9a) 
(9b) 
(9c) 



(10a) 
(10b) 
(10c) 



(The enthalpy change associated with each reaction, AHj, is indicated in brackets). The branching 
of each set of reactions is assumed to be [1:1:1] in these models. Each branch of equations (9) and 
(10) is exothermic, and branches (a) and (b) would require minimal internal re-arrangement. The 
single-product routes in each case are the most exothermic; however, in the case of glycolaldehyde 
formation, the enthalpy change related to reaction branch (10b) is close to that of (10a). While the 
assumed efficiency of branch (c) (i.e parity with the other two branches) may be over-estimated 
in both equations (9) and (10), the production of CO and CH 3 OH (10b) may be very competitive 
with the formation of glycolaldehyde (10a), while it may be less so in the case of methyl formate. 

It thus appears plausible, by this simplistic analysis, that a difference in the efficiencies of 
formation of MF and glycolaldehyde via their primary grain- surface/ice-mantle formation routes 
could account for the order-of-magnitude discrepancy in observational abundances. However, a 
detailed theoretical determination of the branching behavior of these reactions in the solid phase 
would be necessary to confirm this theory. Alternative branching ratios in the photodissociation 
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of methanol, to produce a significant excess of CH 3 over CH 2 OH is also possible, although the 
current - and only - estimate of these values for solid-state methanol, based on experimental rate- 
fitting, suggests precisely the opposite relationship (CH 3 0:CH 2 OH of 1:5; Oberg et al. 2009a). 

A decrease in the efficiency of glycolaldehyde production as described above would be un- 
likely to affect significantly the abundances of other observed species, although the production of 
postulated molecules related to glycolaldehyde, such as dihydroxyacetone (HOCH 2 COCH 2 OH) 
could be affected. The destruction rates of the HCO and CH 2 OH radicals would remain the same, 
while the increased re-formation of CO and CH 3 OH resulting from the alternative product branch 
would have only a minor effect on the already large abundances of those species. 

The glycine formation routes included in this model are sufficiently independent from the 
glycolaldehyde chemistry that glycine production would also be unaffected by the adoption of an 
alternative efficiency for glycolaldehyde formation. Similarly, while acetic acid (CH 3 COOH) is 
a structural isomer of both glycolaldehyde and methyl formate, its different structure means that 
its formation routes are not directly related to those of its isomers; the choice of treatment for 
glycolaldehyde would there have little effect on its own abundance or, through it, that of glycine. 
The chemical behavior of acetic acid is discussed above, in section 3.1.1. 



3.2.3. Formic acid 

Formic acid (HCOOH), shown in Figs. 1-3 panel (b), is formed to a large degree in the 
gas phase, especially for the fast and medium warm-up models; this occurs when formaldehyde 
(H 2 CO) evaporates at ~40 K, reacting with OH to form HCOOH and atomic hydrogen. The 
gas-phase formic acid accretes back onto the grains, and does not evaporate until a temperature of 
around 120 K is reached. A similar or somewhat smaller quantity of the formic acid that ultimately 
evaporates from the grains is formed at low temperature within the ice mantles by the addition of 
HCO and OH, as HCO becomes mobile. 

As such, it is unclear whether the formic acid observed toward hot cores originates in the gas 
phase or on the grains. As discussed above for the case of methyl formate and glycolaldehyde 
production, a plausible alternative product branch for the reaction of HCO and OH exists, yielding 
CO and H 2 (rather than HCOOH); indeed, this branch is somewhat more exothermic than the 
single-product branch. Both branches were included in the model, at a 1 : 1 efficiency ratio. The 
true efficiency for HCOOH formation could therefore be twice as large, or significantly smaller 
than assumed. The rate of the gas-phase reaction (OH + H 2 CO — > HCOOH +H) is reasonably well 
known, albeit only around room temperature (Yetter et al. 1989, and references therein); but the 
amount of formic acid formed as a result is strongly dependent on the amount of formaldehyde 
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that is released into the gas-phase from the grains. 



3.2.4. Formamide 

The behavior of the amino aldehyde formamide (NH 2 CHO), shown in Figs. 1-3 panel (d), 
is somewhat different to that presented in the models of GWH. In this case, the quantity formed 
in the ices, via the addition of HCO and NH 2 radicals, is somewhat larger than previously found. 
However, its peak gas-phase abundance is lower, and falls over time, rather than rising as was found 
by GWH. The latter change in behavior is due to the removal of a single gas-phase reaction from 
the network: NH 2 + H 2 CO — > NH 2 CHO + H. This reaction was present in previous networks with 
a generic collisional rate coefficient of 10" 10 cm 3 s" 1 , but it is likely to have an activation energy 
barrier (based on the barrier of 5.89 kcal/mol (~3000 K) calculated by Li & Lu (2002) for the 
more exothermic product branch NH 3 + HCO). The analogous reaction OH + H 2 CO — >■ HCOOH 
+H, discussed above, has a measured rate of 2 x 10" 12 cm 3 s" 1 , corresponding to a maximum 2% 
efficiency in the branching ratio compared to the alternative H 2 + HCO branch. A theoretical 
study (Li & Lu, 2002) of the NH 2 + H 2 CO ->■ NH 3 + HCO reaction suggested a rate coefficient of 
5.25 x 10" 17 cm 3 s" 1 . For these reasons, the NH 2 + H 2 CO reaction was removed entirely. 

4. Gas-phase routes to glycine formation 

While laboratory experiments have shown success in producing glycine and other amino acids 
by the photolysis of organic ice mixtures (Bernstein et al. 2002; Munoz Caro et al. 2002; Holtom 
et al. 2005; Elsila et al. 2007), gas-phase mechanisms for amino-acid production have also been 
considered. Blagojevic et al. (2003) demonstrated that reaction between acetic acid (CH 3 COOH) 
and ionized or protonated hydroxylamine (NH 23 OH + ) can result in the formation of ionized or pro- 
tonated glycine (NH 23 CH 2 COOH + ) and water. Similar schemes involving propionic (propanoic) 
acid were found to produce the equivalent forms of alanine. The authors estimated rate coefficients 
of > 10" 12 cm 3 s" 1 for each of these processes. 

The most probable source of NH 23 OH + ions in hot cores would be NH 2 OH, either by pho- 
toionization, electron transfer with atomic ions such as He + and C + , or protonation by proton donors 
such as H 3 + , H 3 + or HCO + , all of which processes are already present in the chemical network. 
Protonation by CH 5 + , as suggested by Snow et al. (2007), is not included; however, H 3 + and H 3 + 
are around 2 orders of magnitude more abundant than CH 5 + in the models, while the protonation 
rate coefficients in each case are typically very similar. Protonation rates of hydroxylamine by 
HCO + , H 3 + and H 3 + of (0.59, 1.39 and 6.68) x 10" 9 (r/300 K)" 5 cm 3 s" 1 , respectively, are used 
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in the model. 

Electronic recombination of both ionic and protonated glycine were already added with the 
rest of the glycine-related chemistry; ionization and protonation, with subsequent electronic re- 
combination, are the primary destruction routes for glycine in all the models. Glycine is assumed 
to form in 5 % of recombinations of protonated glycine (in line with values adopted for other 
species and with recent experimental determinations for protonated methanol and dimethyl ether), 
with a generic total recombination rate coefficient of 3 x 10" 7 (r/300 K)" 5 cm 3 s" 1 . The same 
rate is assumed for the recombination of ionized glycine; the ejection of a hydrogen atom is again 
assumed to occur for 5 % of reactions, producing the radical HNCH 2 COOH. This may accrete 
onto the grains to be directly hydrogenated, or may react with a selection of hydrogen-bearing gas- 
phase species, such as HCO, to produce glycine, as outlined in section 2.2 (although the influence 
of these routes on glycine abundance is neglible). Various other gas-phase and grain-surface de- 
struction mechanisms are also available for the HNCH 2 COOH radical, including photodissociation 
and further ion-molecule destruction processes. 

Formation of both hydroxylamine (NH 2 OH) and acetic acid (CH 3 COOH) occurs on the dust- 
grains in this model, and each reaches peak gas-phase abundance at around 130 - 140 K, following 
evaporation. In order to test the gas-phase production of glycine from these species, both of the 
following reactions are added to the network, assuming a rate constant of 10" 11 cm 3 s" 1 for each: 

NH 2 OH + + CH 3 COOH NH 2 CH 2 COOH + + H 2 (11a) 
NH 3 OH + + CH 3 COOH -»■ NH 3 CH 2 COOH + + H 2 (lib) 

Neither of these processes was included in the models already analyzed in section 3. 



4.1. Gas-phase formation results 

With the simple addition of reactions (11), the results are unchanged; the contribution of the 
gas-phase processes is never remotely significant compared to those on the dust-grains, in spite of 
the availability of both reactants at similar times/temperatures. 

It is noteworthy that the abundances of NH 2 OH acheived in the current models are around 3 
orders of magnitude lower than those calculated by GWH; the reason for this is that the major for- 
mation route for NH 2 OH on the grains is the addition of surface/ice-mantle NH 2 and OH radicals. 
The availability of both these reactants is significantly curtailed in the present model by the inclu- 
sion of new hydrogen- abstraction reactions for these species (see Table© and by the more accurate 
consideration of the kinetics of grain-surface reactions that are mediated by activation energy bar- 
riers (see G&P); each reactant is more readily converted back to H 2 or NH 3 before the OH and 
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NH 2 radicals are able to meet. A recent observational search by Pulliam et al. (2012), using the 
NRAO 12 m telescope, detected no NH 2 OH toward a range of previously-observed sources, with 
fractional abundance upper limits of 8 x 10" 12 toward Sgr B2(N). The authors noted that emission 
from a region more compact than 5 arcsec would not be detectable with their instrument, while the 
model suggests that the NH 2 OH would be at least as compact as typical hot-core molecules such 
as methyl formate or ethanol, based on evaporation temperatures. Nevertheless, the new model 
results for NH 2 OH are more in line with the observational data. 

However, a recent laboratory study suggests that the formation of NH 2 OH by successive grain- 
surface hydrogenation of NO is extremely efficient (Congiu et al. 2012). While hydrogenation of 
NO and HNOH are already present in the network, the reaction of H with HNO has until now been 
assumed to lead to H-abstraction, with a barrier of 1500 K. 

To increase the effect of the glycine-forming gas-phase reactions, a parallel branch is added 
to the surface/ice reaction set: H + HNO — > HNOH. A barrier of 1500 K is also assumed; but at 
low temperatures, the reaction competes effectively with hydrogen diffusion, producing a reaction 
efficiency that is close to unity. As a result, large amounts of NH 2 OH are formed, in this case 
during the cold collapse phase, rather than the warm-up phase. A total fractional abundance of 
8 x 10" 7 is present in the ice by the end of collapse, i.e. as the starting point for the warm-up phase 
(see Table|7]). This agrees well with the cold-cloud models presented by Congiu et al. (2012). 

In order to further increase the influence of the gas-phase mechanisms, the rate coefficient of 
reaction (1 la) is increased to 10" 9 cm 3 s" 1 , a plausible maximum value for this reaction. 

The results of this optimized model are shown in Fig. @] The large quantities of NH 2 OH 
present in the dust-grain ices alters the fractional abundances of various species in the model by a 
small degree, although with a negligible direct effect on the production of glycine within the ice. 
The abundance of glycine is in fact slightly reduced, as a result of the somewhat lower quantity 
of acetic acid present in the ice. The destruction of gas-phase acetic acid is rather faster in this 
model, due to reaction (11a). However, when compared to an identical model with reaction (11a) 
switched off (not shown), the contribution of the gas-phase reaction to the total glycine formed in 
this model is around 2 % of the total formed on the grains. Some of the glycine that is formed in 
the gas phase is re-accreted onto the grains, to evaporate later with the rest, while the remainder 
is destroyed by ion-molecule reactions. The immediate impact on gas-phase glycine abundance 
by the newly-included mechanisms may be seen in the small bump at around 130 K; less than 
10" 12 n B is produced. It may be noted that, if no glycine-production mechanisms of any kind were 
active in the ices, and all production thus depended on the gas-phase route, two abundance peaks 
would likely occur, at 130 K and 200 K, corresponding to peak production and peak evaporation, 
respectively. 
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On the the basis that this scenario represents the optimal conditions for gas-phase glycine 
formation to occur in a hot core (including the use of a very optimistic reaction rate), it appears 
that the particular gas-phase formation mechanism tested is unlikely to be significant, especially 
if the mechanisms for formation within the dust-grain ice mantles are in operation as suggested 
by the model. The plausibility of such a large hydroxylamine abundance as would be required for 
significant glycine production is also discussed in Section 5.3. 



5. Spectroscopic modeling 

In order to compare the model results more directly with observations of specific sources, 
a spectroscopic model is constructed that combines observationally-determined temperature and 
density profiles with the temperature-dependent chemical outputs of the model, to simulate emis- 
sion from a selection of key molecules for which observational data exist toward the given source; 
this simulated emission is then convolved to allow direct comparison with the observed spectral 
data. 



5.1. Method 

The warm-up phase of the chemical model produces molecular fractional abundances for pre- 
defined output times, each corresponding to a temperature. The resolution of these output times is 
chosen su ch that the associated te mperature values are separated by less than 1%, over the 8 - 400 



K range. Ivan der Tak et al.l 42000) provide such physical profiles for the envelopes of 14 nearby 
high-mass protostellar sources. Using these temperature-radius and density-radius relations, a 
radius and density is assigned to each temperature value in the chemical model outputs. The 
fractional abundance of each species in the model at each temperature/radius value is multiplied 
by its local density, to produce an absolute abundance value measured in cm" 3 . The resulting 
molecular number-density and temperature profiles are then used to calculate the emission and 
absorption coefficients at each spherically-symmetric radial position in the core. 

Partition functions and emission-line data from the JPL line list and the Cologne Database 
for Molecular Spectroscopy are imported into the code as needed, using the 'Splatalogue' online 
database^ The spectral coefficients may be calculated by this method to incorporate all molecular 
emission at any given frequency, allowing line blends and line confusion to be taken into account. 
Gaussian line profiles are assumed, using observationally-determined line widths. A channel width 
(frequency-bin size) much smaller than the line width is chosen (in general, 0.1 MHz or ~0.08 
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km/s), to ensure that the lines are well resolved. 

Under conditions of local thermodynamic equilibrium, the equation of radiative transfer may 
be parameterized thus (Rohlfs & Wilson 2006): 

I u (s) = I u (0) e~ T ' As) + [ A B v (T(r)) e~ T " dr (12) 
Jo 

where l v is the specific intensity, r v is the optical depth at frequency v, T is temperature, s is the 
line-of-sight distance, and B U (T) is the Planck function (equal to ejn v , under LTE). Using the grid- 
ded, spherically- symmetric input data, this equation is integrated numerically along lines of sight 
(using Simpson's rule), beginning directly on-source, and moving outward through parallel lines 
of sight until a user-defined maximum offset is reached. Line-of-sight integrals of the absorption 
coefficient similarly yield the optical depths. This procedure results in an intensity map of the hot 
core for each frequency channel. 

The final step is to convolve the emission with a Gaussian beam centered on the source, using 
a beam-size appropriate to a chosen instrument observing at a given frequency. The intensity map 
is numerically integrated with the Gaussian function, over the radius of the hot-core envelope, out 
to a few beam widths. The technique automatically accounts for beam-dilution, and makes use 
of the precise emission structure to calculate an intensity value at each frequency. The resulting 
simulated spectrum may be directly compared with main-beam temperature (r mb ) spectra obtained 
from the telescope for which the convolution has been produced. 



5.2. Spectroscopic simulations of NGC 6334 IRS1 



The molecular mm/sub-mm line emissio n from a selection of the sources observed by Ivan der Tak et al 



(120001) was surveyed by lBisshop et al.l(|2007l) using the James Clerk Maxwell and IRAM 30 m tele- 
scopes, providing information on dozens on complex organic species. A member of this subset, 
the hot-core source NGC 6334 IRS1, is relatively closeby (1.7 kpc), and exhibits relatively narrow 
(AV = 5 km/s), strong emission lines (see van der Tak et al. 2000 and Bisschop et al. 2007, and 
references therein). It is also well situated in the sky (Dec. (1950) = -35°44') for future observa- 
tion with the ALMA telescope (latitude -23°). For these reasons, it may be identified as a plausible 
candidate for the future detection of glycine. The spectroscopic model described above is therefore 
used to assess the potential for the detection of glycine in this source, on the assumption that its for- 
mation and subsequent behavior are well described by the chemical model. The physical profiles 



2 splatalogue.net 
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provided by van der Tak et al. (2000) are used. In the case of glycine, only emission lines whose 
frequencies are within the observable ranges of the JCMT, IRAM 30 m and ALMA telescopes are 
modeled. 

Fig. \5\ shows - for illustrative purposes - the predicted emission from this source over the 
241.350 - 241.425 GHz range, using a channel spacing of 0.488 MHz. These spectral calculations 
consider emission from all species in the model for which spectroscopic data exist, using only 
ground-state transitions. Glycine emission, at 241.373 GHz, is highlighted in red; glycine con- 
former I is simulated, but not conformer II, due to uncertainty in the partition functions. The upper 
panel of the figure shows the unconvolved emission for an on-source line of sight - the spectrum 
may be considered as that expected from a pencil beam of infinite resolution directed at the source. 
Individual lines are identified in the upper panel; an identification is made where more than half of 
the emission within a channel derives from a single molecule. The emission from glycine is seen 
to be slightly blended with emission from glycolaldehyde (CH 2 (OH)CHO) in the line wings, but 
nevertheless shows a clear peak. The middle panel shows the spectrum after convolution of the 
emission from the entire source (same temperature scale) using a beamwidth of 0.4 arcsec, which 
encompasses the strongest predicted emission region for glycine (see below) and which should be 
achievable with ALMA at these frequencies (although the simulation is for a single-dish telescope, 
rather than an array). The glycine peak is preserved, with a line strength on the order of 1 K. 
The lower panel of the figure shows the predicted spectrum (different temperature scale) using a 
beamwidth appropriate to the James Clerk Maxwell Telescope at these frequencies (~10.4 arcsec); 
the glycine emission line strength is on the order of 1 mK. 

A selection of the strongest emission lines from methanol, methyl formate and glycine have 
been simulated using the method above. Methanol and methyl formate were both detected by 
Bisschop et al. (2007) toward this source; both are typically-observed hot-core species. In these 
calculations, unlike those illustrated in Fig. [51 only the emission from an individual molecule is 
considered at the chosen line fr equency; this e nable s a more straight-forward comparison with 



the unblended lines observed by iBisshop et al.l (|2007|) . and makes the calculations less computa 



tionally intensive. Furthermore, while line identifications may be made fairly easily in the uncon- 
volved "pencil-beam" simulations, assessing the contributions made by each molecular line to the 
convolved emission spectra is much more time-consuming. The consideration of the detectability 
of each glycine line in the following analysis therefore does not consider potential line blends - 
simply the expected line strength. More detailed assessments of the expected line emission from 
various species will be made in future work. 

Fig. [6] shows the emission profiles predicted for the line center of a strong methyl formate 
line (left panel) and a strong glycine line of similar frequency (right panel; glycine emission is 
also shown in left panel for direct comparison). The local emission strength (i.e. assuming a 
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pencil beam) is shown in terms of brightness temperature for lines of sight targeted on radial 
displacements from the source position on the sky. Figure [7J shows the same profiles mapped 
onto sky coordinates; these emission maps are then convolved with the beam to produce simulated 
detected brightness temperatures. The most significant emission from methyl formate is seen to 
emanate from a region of diameter ~7000 AU, or 4 arcsec. All potentially-detectable emission 
from glycine comes from a region ~0.8 arcsec in diameter, although the most significant emission 
region is of diameter ~0.4 arcsec, corresponding to a radius of ~1000 AU. 



5.2.1. Comparison with observations 

Tables |9l [101 and [TT1 show details of the predicted emission from each molecule, using model 
outputs corresponding to the fast, medium and slow warm-up models, respectively. Line frequen- 
cies are shown in GHz, along with upper-level energies in K; in some cases, each 'line' represents 
two or more lines of identical or very similar frequencies and/or spectral characteristics. The 
column marked 'peak simulated local intensity' indicates the maximum local emission strength 
(assuming a pencil beam). Beam widths represent the full-width half-maximum of the Gaussian 
beam used; the precise beam width appropriate to the JCMT or IRAM 30 m is calculated by in- 
terpolation of telescope- specific data available online. The convolved intensity of the line is the 
relevant quantity for discerning the detectability of a line; however, for comparison with the ob- 
served lines, the integrated intensity (in units of K km s" 1 ) provides a better test, as this quantity 
removes the dependence on the average line width (AV); this allows all of the emission predicted 
by the model to be compared against all of the emission detected in the observations. The final col- 
umn of these tables shows the integrated intensities measured for the methanol and methyl formate 
lines by Bisschop et al. (2007). 

It may be seen immediately that the simulated integrated intensities of the methanol lines are 
an excellent match to the observed values, being generally around a factor of 2 lower, and showing 
only minor variation across the different warm-up timescale models. The simulations also provide 
a good match regardless of the upper-level energies and beam sizes (which vary according to 
frequency). The methyl formate emission is much more variable between model timescales; the 
fast model results provide an excellent match to all integrated intensities, falling within a factor 
less than 1.5 of the observed values. The medium timescale model produces values that are around 
2-3 times lower than those observed, similar to the agreement shown with the methanol lines. 
However, the slow warm-up timescale results produce integrated intensities only one hundredth 
part as large as presented in the observations of Bisschop et al. (2007). In this case, methyl 
formate suffers rapid gas-phase destruction, as described in Section 3.2.2, producing much lower 
instantaneous fractional abundances and much faster decline over time/temperature. The two faster 
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warm-up models therefore appear the most favorable in comparison with observations. 



5.2.2. Glycine predictions 

The tables show predictions for six of the strongest glycine lines, or suites of lines. Each 
has an upper-level energy sufficiently low to be well-populated at the 200+ K temperatures at 
which glycine fractional abundances are strongest in the models. The predicted intensities between 
different lines (within each model) are very consistent, and peak local intensities range from around 
1 K to 10 K, corresponding to fast and slow warm-up timescales, respectively. It is apparent that 
even for this very nearby and otherwise favorable hot-core source, the glycine emission strength 
is predicted to be very moderate; convolving the emission with a beam appropriate to the JCMT 
(~20 - 22 arcesec) produces line strengths of no more than 10 mK in the most optimistic model. 
In the case of the fast warm-up model, which reproduces methyl formate emission most accurately, 
the glycine line strengths are no more than 250 fxK. 

Convolution of the emission using a beam appropriate to the IRAM 30 m telescope for each 
line (~10 - 12 arcsec) produces line strengths of around 1, 3, and 30 mK, for fast, medium, and 
slow models respectively. This suggests that one of the strongest glycine emission lines could 
plausibly be detected toward the given source using this instrument, assuming that the slow model 
values are accurate; a 3-a line detection would require a noise level of around 10 mK (assuming no 
line blending, which is not modeled here). However, detection assuming either of the other model 
values would require impractical signal-to-noise values. 

In order to make an approximate estimate of the detectability of glycine using ALMA, the 
same single-beam convolution method is used, along with two plausible, fixed beam sizes: a 
generic 1 arcsec beam, and a 0.4 arcsec beam to match the size of the region of strongest glycine 
emission (see Fig. [6]). The reduced beam-dilution produced in these cases results in much greater 
line strengths, and suggests that the detection of a strong glycine line with ALMA is very plausible 
(line blending and model inaccuracies notwithstanding). The weakest emission predicted for any 
of the simulated glycine lines, with a beam width of 1 arcsec, is ~87 mK, while at 0.4 arcsec the 
weakest line has a strength of ~350 mK. Under the optimistic slow warm-up model conditions, 
convolved line strengths are close to 10 K for the strongest lines, using the smaller beam. For 
comparison, the ALMA on-line sensitivity calculated indicates that, with 32 antennae (Cycle 1) 
and a resolution of 0.4 arcsec, ALMA would achieve an rms of 200 mK per 1 km/s velocity bin in 
1 hour. 



http://almascience.eso.org/call-for-proposals/sensitivity-calculator 
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The line simulations therefore suggest that while detection of glycine using single-dish in- 
struments is unlikely, even for this very favorable source, detection with ALMA is highly plausible 
(ignoring possible line-blending effects). 



5.3. NH 2 OH toward W3 IRS5 

Pulliam et al. (2012) recently looked for hydroxylamine emission toward a selection of 
sources including Sgr B2(N), Orion KL and W3 IRS5. They found no NH 2 OH emission using 
data obtained with the NRAO 12 m telescope, with noise levels of a few mK. Both the GWH 
model and the enhanced NH 2 OH model of section 4.1 exhibit large fractional abundances for this 
molecule, albeit with different formation mechanisms in each case. In order to test whether the 
model values are consistent with the non-detections of Pulliam et al., the spectroscopic model is 
again applied to the chemical model results. Out of the sources observed by both Pulliam et al. 
(2012) and van der Tak et al. (2000), the only source present in both datasets is W3 IRS5. For this 
simulation, a generic line width of 9 km/s is assumed, following Pulliam et al. (2012). 

Tables [T2l and [T3l show the results of a selection of NH 2 OH emission-line simulations for W3 
IRS5, using, respectively, the standard medium warm-up timescale model results, and the results 
of the model of Sec. 4.1 that employs a gas-phase glycine formation mechanism and a large initial 
NH 2 OH abundance on the grains, produced by the low-temperature formation of NH 2 OH on dust 
grains during the collapse phase. The results shown in Table [12] are all consistent with the non- 
detection; Pulliam et al. (2012) find a l-a noise level of 4.3 mK in their W3 IRS5 data, while 
the simulated intensities are of order 10 /iK or less. However, the model with enhanced NH 2 OH 
formation produces maximum convolved line intensities of around 50 mK, or around 1 order of 
magnitude greater than the noise level. 

Since the simulated line strengths scale well with the peak fractional abundance of NH 2 OH 
obtained in the chemical models, the observed values would require a peak gas-phase abundance of 
< 10" 7 for the models to reproduce the non-detection of hydroxylamine, assuming all other param- 
eters in the models are correct. The implication for the models is that the cold-stage conversion of 
NO/HNO to NH 2 OH on the grains is too efficient, nominally by a factor of 10. This would suggest 
that the gas-phase formation of significant quantities of glycine is also less probable, in all sources. 
The apparent over-estimate of hydroxylamine production could be explained by several possibil- 
ities; one is that the chemical network does not include a sufficent number of competing reaction 
mechanisms for all of the states of hydrogenation between NO and NH 2 OH. Another possibility is 
that the determination that the reaction H + HNO — > H 2 NO has an insignificant barrier, by Congiu 
et al. (2012), is dependent on experimental conditions that are inappropriate to these models, such 
as the ice composition, which is predominantly water in the ISM. Alternatively, the hydroxylamine 
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abundance in the ice may be adversely affected either by photo-destruction or hydrogen abstraction 
processes within the ice, whose rates or barriers may be inaccurately quantified in the models. 



6. Discussion 

6.1. Dust-grain chemistry 

The general behavior of complex molecules as simulated in the new models is broadly sim- 
ilar to that seen in the previous models, such as G&H and GWH; evaporation of grain-surface 
molecules into the gas-phase occurs around the same temperatures and in similar quantities. How- 
ever, the formation temperatures of molecules on the grains are somewhat different, due both to 
the distinction between surface and ice-mantle populations and to the difference in diffusion rates 
in each medium; HCO-related species are formed strongly around 25 K, while CH 3 -related species 
are formed around 20 K - both som ewhat lower than p revious models, though more in line with 
recent experimental results such as lobergetal.l(b009al) . 

The interaction between the surface and ice-mantle chemistry in the models is complex; mo- 
bile radicals may land on the ice surface from the gas phase or may diffuse out of the mantle, 
to react to form more complex species, which show a net diffusion back into the mantle as more 
volatile species move out. The migration of mantle species to the surface as temperatures increase 
tends to mitigate the 'trapping' effect that would be expected using a model in which there were 
no surface-mantle diffusion. Much complex chemistry also occurs within the mantles themselves, 
particularly at elevated temperatures. The importance of mantle versus surface processes is de- 
pendent on the source of the reactants, the mobilities of the reactants, the temperature, and the 
general composition of the surrounding ice; no general statement can therefore be made on this 
point other than that the diffusion within the ice mantles, and between the surface and the bulk 
ice, appears to play a major part in hot-core chemistry. (Although the importance for specific, 
individual molecules is discussed in section 3). 

Aside from the effects of the differing treatment of the physical structure of the ices employed 
in this model, the inclusion of a more comprehensive set of hydrogen-abstraction reactions in the 
ice chemistry also demonstrates the key role of water and OH in producing other free radicals 
within the ice that may go on to form more complex organic structures. The OH radical is highly 
reactive, and may abstract hydrogen from stable molecules with only a small activation energy - 
typically lower than the barrier to OH diffusion. Thus, hydrogen abstraction by OH may occur 
with a net rate equal to its rate of diffusion, in the case where the target molecule is less mobile 
than OH. The large quantities of water present in the ice makes OH the most commonly produced 
photodissociation product, caused by the weak, cosmic ray-induced, UV radiation that permeates 
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the hot core. 

The complex molecules, such as glycolaldehyde (CH 2 (OH)CHO), ethanol (C 2 H 5 OH) and 
methyl formate (HCOOCH 3 ), formed by the addition of HCO or CH 3 to the photoproducts of 
methanol (CH 3 OH) are produced early, at low temperatures, before OH becomes mobile (or abun- 
dant). Their production is therefore reflective of the direct photodissociation of methanol into 
CH 2 OH, CH3O and CH 3 . 

While greater temperatures increase OH mobility, below 40 K its primary destruction route 
(and the primary formation route for HCO) is the reaction OH + H 2 CO — > HCO + H 2 0, which 
is assumed to require no activation energy. The mobility of H 2 CO itself drives the reaction, and 
it is not until H 2 CO evaporates at around 40 K that OH becomes important in providing radicals 
other than HCO. After this, HCO becomes less prevalent in the ices, while CH 3 is still produced 
in some abundance by methanol photodissociation, resulting in a moderate bias toward methyl- 
group addition over aldehyde production. At this point, ammonia (NH 3 ) and methanol, as major 
constituents of the ice, become the dominant reaction partners for OH, producing NH 2 and CH 2 OH 
(CH 3 is also produced, but the higher activation energy for this process renders it uncompetitive 
versus the hydroxymethyl channel). 

It is therefore at temperatures above ~40 K that the formation of amino-organics becomes 
most active, by addition of NH 2 to methanol products, or to more complex radicals. However, 
at even higher temperatures, NH 2 also becomes an important source of radicals in itself, as the 
barriers to hydrogen abstraction by NH 2 are relatively low compared to its diffusion barrier, as is 
the case for OH. Above around 55 K, the photodissociation of water into OH typically results in 
the subsequent conversion of ammonia to NH 2 , which becomes the ultimate agent of hydrogen 
abstraction from organic molecules in the ice. This continues until ammonia, the source of NH 2 
radicals, evaporates at ~120 K. 



6.2. Glycine formation 

The inclusion in the model of ice surface/mantle formation mechanisms for glycine indicates 
that it may be formed in a number of ways, whose degree of influence depends on the timescale of 
the hot-core evolution. The abundance of acetic acid (CH 3 COOH) is perhaps most important in de- 
termining which mechanism is most important. In the slow warm-up model, acetic acid abundance 
is relatively high within the ices, and the abstraction of hydrogen by NH 2 , followed by the addition 
of another amine group, is the most important glycine-formation route. The same is true, albeit 
to a lesser extent, for the medium warm-up timescale model. The fast warm-up model produces 
less acetic acid, and the addition of CH 2 NH 2 and HOCO - products of OH-induced hydrogen ab- 
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straction from methylamine and formic acid, respectively - becomes the major formation route 
in conjunction with OH addition to the NH 2 CH 2 CO radical, which is formed by abstraction from 
the amino aldehyde glycinal (NH 2 CH 2 CHO). It is questionable as to whether the high acetic acid 
abundances produced by the slow model are supported by observations, indicating both that glycine 
abundances may be lower than suggested by that model, and that no single ice surface/mantle pro- 
cess of the three that are tested in this work should be strongly dominant. 

It may also be noted that, in the case of the CH 2 NH 2 + HOCO glycine-formation route, a 
plausible alternative channel, producing CH 3 NH 2 + C0 2 , may also be important. In this model a 
branching ratio of 1: 1 is assumed; however, the accuracy of this assumption cannot be ascertained 
without experimental or computational data for the solid-phase process (as opposed to the gas- 
phase equivalent, in which glycine production would be prohibitively disfavored). The effect of a 
lower efficiency for glycine production through this method would be strongest in the fast warm-up 
scenario, although the availability of alternative routes means that glycine abundance is unlikely 
to fall by more than a factor of 2 as a result; conversely, a more favorable branching ratio could 
increase the glycine production by a similar factor under fast warm-up conditions. 

A gas-phase mechanism for the formation of glycine was also tested in this model, consisting 
of the addition of protonated hydroxylamine (NH 3 OH + ) and acetic acid (CH 3 COOH), followed 
by electronic recombination. Even assuming the highest plausible rate, with large quantities of 
hydroxylamine and acetic acid produced on the grains, the amount of glycine produced is only on 
the order of 10" 12 n H > or 1% of the amount formed on the dust grains. The limited comparison of 
simulated NH 2 OH emission with observations suggests that the actual fractional abundance may 
indeed be lower than assumed in that model. It is possible that the fraction of dissociative recom- 
bination of protonated glycine that produce glycine itself may be somewhat higher than assumed 
(5%), but it is unlikely to be underestimated by more than a factor of a few. The combination of 
favorable conditions and rates that would be required therefore suggest that gas-phase production 
by this method is unlikely to produce significant quantities of glycine. 

The possibility, demonstrated in the models, for significant glycine production to occur through 
several distinct chemical routes suggests that the formation of glycine within the ice mantles is no 
less plausible than the formation of more commonly observed complex organic molecules. 



6.3. Detectibility of glycine 

The peak gas-phase fractional abundances for glycine of 8 x 10" 11 - 8 x 10" 9 produced by the 
chemical models suggest that this molecule could be reasonably abundant in regions within hot 
cores. However, the LTE radiative transfer and convolution models suggest that the detection of 
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glycine toward a favorable source may be impossible with single-dish instruments, assuming any 
but the most optimistic of the chemical models; this is due largely to the extremely compact emis- 
sion exhibited by glycine. The model findings are therefore in agreement with the current lack of a 
detection of this molecule. However, it appears that the new ALMA telescope will have sufficient 
sensitivity and spatial resolution to easily detect the strongest glycine emission lines, assuming 
that the emission is not obscured by that from other molecules, and that the source is sufficiently 
closeby that the highly compact glycine-emission region is well resolved. Detection will likely 
require the careful choice of sources that are both nearby and that exhibit narrow emission lines. 
For this reason, detection toward the well-observed Galactic Center source Sgr B2(N) may be 
questionable. While ALMA should allow the distinction of different structures within that source, 
which contribute to the considerable line widths (Belloche et al. 2008), it is not clear whether a 
predicted compact glycine-emission region (of temperature ~200 K) could be resolved sufficiently 
well both to produce high enough line strengths for detection and to filter out contamination from 
other molecules. No simulations - such as those produced above for other sources - are possible, 
due to the absence of both a temperature and density profile for Sgr B2(N) in the literature. 

The most important factor in the detectibility of glycine toward hot cores is likely to be its 
binding energy, which determines its evaporation temperature, and thus its radius of peak fractional 
abundance. The simulations of glycine emission in NGC 6334 IRS1 suggest that the strongest 
glycine emission should emanate from a region ~1000 AU in diameter within that source. This 
putative emission is predicted to be detectable with ALMA. The single-dish convolution of the 
simulated emission shown in Section 5.2, using a beam width of 0.4 arcsec as appropriate to 
ALMA, represents a first approximation to the detectibility of glycine at that resolution; the use 
of an interferometer would produce spatial filtering such that the more extended emission of other 
hot-core molecules could be resolved out, lessening the impact of contamination from other lines 
on the glycine emission. 



6.4. Combined chemical/spectroscopic models 

The combined chemical and spectroscopic models produce a very good match with observed 
methanol line emission toward NGC 6334 IRS1. The match for methyl formate is also good, in the 
case of the two shorter-timescale models; warm-up timescales of around 1 Myr may not therefore 
be appropriate to hot cores, although the precise binding energy of methyl formate is important 
to this determination. The methyl formate results also indicate that, while other mechanisms may 
plausibly exist for the formation of methyl formate in hot cores, the addition of methanol photo- 
fragments within the dust-grain ice mantles, as modeled here, appears adequate to explain observed 
abundances. Meanwhile, the structural isomer of methyl formate, glycolaldehyde, is almost cer- 
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tainly over-produced, although this may be remedied by a more careful consideration of the product 
branching ratios of the reaction between HCO and CH 2 OH radicals in the ices. Branching ratios 
favoring an alternative [CO + CH 3 OH] channel would seem plausible. 

All of the radiative transfer calculations presented here are highly dependent on the observationally- 
determined temperature-density-radius profiles provided by van der Tak et al. (2000). While the 
simulations of methanol and methyl formate emission agree remarkably well with recent observa- 
tional data, the uncertainty in the profiles naturally increases at smaller radii, and the fidelity of 
a spherically-symmetric model must inevitably break down where outflows or clumpy structure 
become important. Nevertheless, in the absence of more precise data, perhaps also to be provided 
by ALMA, the simple models used here provide a first quantitative estimate of the behavior and 
detectability of glycine in star-forming regions. A more comprehensive assessment would also in- 
corporate the analysis of potential line blending, using spectroscopic models that include all species 
for which data exist, as well as an explicit treatment to simulate interferometric observations where 
appropriate. 

7. Conclusions 

The main conclusions of this work are summarized below: 

1. The hot-core chemical models predict peak gas-phase glycine abundances of ~8 x 10" 11 - 
8 x 10" 9 nu (dependent on warm-up timescale). The peak abundance values are attained at 
around 200 K, when glycine evaporates from the dust grains. 

2. Glycine is found to form almost exclusively within or upon dust-grain ice mantles, beginning 
at temperatures of ~40-50 K, and ending at ~ 120 K. The three main radical-radical addition 
mechanisms investigated appear to have approximately equal influence on the resultant gas- 
phase quantities. 

3. Gas-phase glycine formation involving acetic acid (CH 3 COOH) and protonated hydroxy - 
lamine (NH 3 OH + ) is insignificant under even the most optimistic conditions. 

4. Related organic species such as propionic acid, propionaldehyde (propanal), and the amino 
aldehyde glycinal are predicted also to attain significant gas-phase abundances, although 
their detectibility would be dependent also on their emission radii, via their evaporation 
temperatures. 

5. The production of glycinal (NH 2 CH 2 CHO) and acetic acid (CH 3 COOH) - two key molecules 
in the formation of glycine in the ice mantles - are strongly dependent on the earlier produc- 



-35- 



tion of acetaldehyde (CH 3 CHO). Acetaldehyde and acetic acid could provide observational 
information on potential glycine abundances in the absence of direct detection. 

6. The current lack of a detection of glycine toward hot cores may be explained by its high 
temperature of evaporation, which results in a small emission radius. Glycine is not expected 
to be detectable with single-dish instruments over realizable integration times toward any 
hot-core source. 

7. The spectroscopic model suggests glycine may be detectable using ALMA at sub-arcsecond 
resolutions toward bright, nearby sources with relatively narrow emission lines, such as NGC 
6334 IRS 1. 

8. Observed methyl formate abundances are well reproduced by the two shorter-timescale mod- 
els, indicating that grain-related formation mechanisms are sufficient (as suggested by previ- 
ous models), and that (8-400 K) warm-up timescales of around 1 Myr may be inappropriate 
for hot cores. 

9. The anomalously high abundance produced by the models for the structural isomer of methyl 
formate, glyclolaldehyde, may be remedied by the consideration of alternative branches for 
the HCO + CH 2 OH reaction in the ice mantles, favoring the formation of CO + CH 3 OH. 

10. Cosmic ray-induced photodissociation of water ice, to produce OH, is crucial to the produc- 
tion of radicals within the bulk ices. The OH radicals easily abstract hydrogen from other 
molecules, including ammonia, to produce reactive molecular radicals that may combine to 
form more complex species. 

1 1 . NH 2 , formed by hydrogen abstraction from ammonia by OH, also appears to be an important 
agent of radical propagation within the ice mantles. This process, and the related production 
of amino-organics, becomes most significant at temperatures between ~40-120 K. The ob- 
served nitrogen content of hot-core molecules may thus be related to the timescale spent in 
this temperature range. 

12. Much, if not most, of the chemistry that produces complex organic molecules occurs within 
the ice mantles, while the more rapid diffusion that obtains on the ice surface allows different 
reaction kinetics to occur concurrently. Thus, the explicit treatment of these two phases, their 
coupling via thermal diffusion, and the interaction of the dust/ice surface with the gas phase 
are all essential elements in the chemical model. 

13. The simple spectroscopic model used here provides a direct comparison between chemical 
models and observational spectroscopic data for complex molecules in specific sources, an 
element that has been lacking until now in the study of hot cores. 
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Table 1 . Initial fractional abundances, with respect to total hydrogen, of elements and H 2 at the 

start of the collapse phase. 



Species, i 


Ui/riH t 


H 2 


0.499 


H 


2.0(-3) 


He 


0.09 


C 


1.4(-4) 


N 


7.5(-5) 





3.2(-4) 


s 


8.0(-8) 


Na 


2.0(-8) 


Mg 


7.0(-9) 


Si 


8.0(-9) 


P 


3.0(-9) 


CI 


4.0(-9) 


Fe 


3.0(-9) 



U(B) = A 
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Table 2. A selection of barrierless surface/mantle radical-addition reactions relevant to glycine 
or related species included in the chemical network 



Reaction 



HOCO 


+ 


CH2NH2 


-> 


NH 2 CH 2 COOH 












-> 


CH3NH2 


+ 


C0 2 


NH 2 


+ 


CH2COOH 


->• 


NH2CH2COOH 






OH 


+ 


NH 2 CH 2 CO 


-> 


NH 2 CH 2 COOH 






H 


+ 


HNCH 2 COOH 


-> 


NH 2 CH 2 COOH 






CH 2 OH 


+ 


HNCH 2 COOH 




NH 2 CH 2 COOH 


+ 


H 2 CO 


CH 3 


+ 


HNCH 2 COOH 


-> 


NH 2 CH 2 COOH 


+ 


H 2 CO 


HCO 


+ 


HNCH2COOH 




NH2CH2COOH 


+ 


CO 


HOCO 


+ 


HNCH 2 COOH 




NH 2 CH 2 COOH 


+ 


C0 2 


NH 


+ 


CH2COOH 


-> 


HNCH2COOH 






HCO 


+ 


CH 2 NH 2 


— > 


NH 2 CH 2 CHO 












— >■ 


CH3NH2 


+ 


CO 


NH 2 


+ 


CH 2 CHO 


— > 


NH 2 CH 2 CHO 






H 


+ 


NH 2 CH 2 CO 




NH 2 CH 2 CHO 






CH 2 OH 


+ 


NH 2 CH 2 CO 


— > 


NH 2 CH 2 CHO 


+ 


H 2 CO 


CH3O 


+ 


NH 2 CH 2 CO 


— >■ 


NH 2 CH 2 CHO 


+ 


H 2 CO 


HCO 


+ 


NH 2 CH 2 CO 




NH 2 CH 2 CHO 


+ 


CO 


HOCO 


+ 


NH 2 CH 2 CO 




NH 2 CH 2 CHO 


+ 


C0 2 


H 


+ 


HNCH 2 CHO 


-> 


NH 2 CH 2 CHO 






CH 2 OH 


+ 


HNCH2CHO 




NH 2 CH 2 CHO 


+ 


H 2 CO 


CH3O 


+ 


HNCH2CHO 




NH 2 CH 2 CHO 


+ 


H 2 CO 


HCO 


+ 


HNCH 2 CHO 


— ^ 


NH 2 CH 2 CHO 


+ 


CO 


HOCO 


+ 


HNCH2CHO 


— >■ 


NH 2 CH 2 CHO 


+ 


C0 2 


NH 


+ 


CH 2 CHO 




HNCH 2 CHO 






CH3CO 


+ 


OH 


— > 


CH3COOH 






CH 3 


+ 


HOCO 


— >■ 


CH3COOH 






H 


+ 


CH 2 COOH 


— >■ 


CH3COOH 






CH 2 OH 


+ 


CH2COOH 


— > 


CH3COOH 


+ 


H 2 CO 


CH3O 


+ 


CH2COOH 




CH3COOH 


+ 


H 2 CO 


HCO 


+ 


CH 2 COOH 


— t 


CH3COOH 


+ 


CO 


HOCO 


+ 


CH2COOH 




CH3COOH 


+ 


co 2 


CH 2 


+ 


HOCO 




CH 2 COOH 






CH 2 


+ 


CH 2 COOH 




C2H4COOH 






CH 3 


+ 


CH2COOH 


— > 


C2H5COOH 






H 


+ 


C2H4COOH 


— ► 


C2H5COOH 






CH 2 OH 


+ 


C2H4COOH 


— > 


C2H5COOH 


+ 


H 2 CO 


CH3O 


+ 


C2H4COOH 




C2H5COOH 


+ 


H 2 CO 


HCO 


+ 


C2H4COOH 


— > 


C2H5COOH 


+ 


CO 


HOCO 


+ 


C2H4COOH 




C2H5COOH 


+ 


co 2 


HOCO 


+ 


C2H5 


— > 


C2H5COOH 












— > 


C2H6 


+ 


co 2 


OH 


+ 


C2H5CO 


— >■ 


C2H5COOH 






CH 3 


+ 


CH 2 CHO 


— »• 


C2H5CHO 






H 


+ 


C2H4CHO 


->• 


C 2 H 5 CHO 






CH 2 OH 


+ 


C2H4CHO 




C2H5CHO 


+ 


H 2 CO 


CH3O 


+ 


C2H4CHO 




C 2 H 5 CHO 


+ 


H 2 CO 
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Table 2 — Continued 



Reaction 


HCO 


+ 


C2H4CHO - 


-> C2H5CHO 


+ 


CO 


HOCO 


+ 


C,H 4 CHO - 


■> CjHsCHO 


+ 


C0 2 


CH 2 


+ 


CH 2 CHO 


4 C2H4CHO 






HCO 


+ 


C 2 H 5 


-> C2H5CHO 






H 


+ 


C 2 H 5 CO 


-> C2H5CHO 






CH 2 OH 


+ 


C 2 H 5 CO 


-> C2H5CHO 


+ 


H 2 co 


CH3O 


+ 


C 2 H 5 CO 


C2H5CHO 


+ 


H 2 CO 


HCO 


+ 


C 2 H 5 CO 


-> C2H5CHO 


+ 


CO 


HOCO 


+ 


C 2 H 5 CO 


-> C2H5CHO 


+ 


C0 2 



Note. — Not all reactions involving each species are shown. Equal 
branching ratios are assumed, where multiple channels are given. The 
formation of the relevant radicals may occur through photodissocia- 
tion or through hydrogen-abstraction reactions, detailed in Table|5] 
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Table 3. Activation energy barrier information for important barrier- mediated surface/mantle 

reactions 



# 


Reaction 


Ea 


Width 


Ref. 






(K) 


(A) 




1 


H + CO — y HCO 






Andersson (priv. comm.) 


2 


H + H2CO — y CH2OH 


4500 


1.35 


Goumans (priv. comm.) 


3 


H + H 2 CO -S- CH3O 


2320 


1.35 


Goumans (priv. comm.) 


4 


H + H2CO — y HCO + H2 


2960 


1.22 


Goumans (priv. comm.) 


5 


H + CH3OH — ¥ CH 2 OH + H 2 


4380 


1.0 


Goumans (priv. comm.) 


6 


H + CH3OH — y CH3O + H 2 


6640 


1.0 


Goumans (priv. comm.) 


7 


H + CH4 — y CH 3 + H 2 


5940 


2.17 


Andersson et al. (2009); Garrod & Pauly (2011) 


8 


H + CH3CHO — s- H 2 + CH3CO 


2120 


1.0 


Warnatz (1984) 


9 


H + HCOOCH3 — ¥ H2 + CH3OCO 


3970 


1.0 


Good & Francisco (2002) 


10 


H + HCOOCH3 — > H2 + HCOOCH9 


4470 


1.0 


Good & Francisco (2002) 


11 


H + HCOOH — ¥ H 2 + COOH 


6150 


1.0 


Evans-Polanyi, #8/#12 


12 


H + NH2CHO —y H2 + NH2CO 


3500 


1.0 


Syrstad & Turecek (2001) 


13 


H + CH 2 (OH)CHO — y H 2 + CH 2 (OH)CO 


2350 


1.0 


Evans-Polanyi, #8/#12 


14 


H + HNCO — y NH 2 + CO 


2300 


1.0 


Tsang (1992) 


15 


H + HNCO — y NH2CO 


1390 


1.0 


Nguyen et al. (1996) 


16 


H 2 + OH — y H 2 + H 


2100 


1.0 


Baulch et al. (1984) 


17 


CH2OH + CH3CHO — y CH3OH + CH3CO 


5300 


1.0 


Evans-Polanyi 


18 


CH 2 OH + CO — y CH 2 (OH)CO 


2500 


1.0 


Garrod et al. (2008) 


19 


CH2OH + H2CO — y CH3OH + HCO 


2940 


1.0 


Tsang H9871 


20 


CH2OH + NH2CHO — y CH3OH + NH2CO 


4320 


1.0 


Evans-Polanyi 


21 


CH 2 OH + CH 2 (OH)CHO — y CH 3 OH + CH 2 (OH)CO 


3100 


1.0 


Evans-Polanyi 


22 


CH2OH + HNCO — y CH3OH + OCN 


1200 


1.0 


Estimate 


23 


CH2OH + C2H5CHO — y CH3OH + C2H5CO 


3100 


1.0 


(#21) 


24 


CH 2 OH + C3H4 — y CH3OH + C3H3 


4960 


1.0 


Evans-Polanyi 


25 


CH2OH + C3H6 — y CH3OH + C3H5 


4790 


1.0 


Evans-Polanyi 


26 


CH 2 OH + CH 2 NH -y CH3OH + H 2 CN 


1200 


1.0 


(#22) 


27 


CH 2 OH + CH3NH2 —y CH3OH + CH 2 NH 2 


6780 


1.0 


Evans-Polanyi 


28 


CH2OH + CH3CN — y CH3OH + CH2CN 


6200 


1.0 


Evans-Polanyi 


29 


CH 2 OH + C 2 H 5 CN — y CH3OH + C 2 H 4 CN 


6100 


1.0 


Evans-Polanyi 


30 


CH 2 OH + C3H7CN ->■ CH3OH + C 3 H 6 CN 


6100 


1.0 


(#29) 


31 


CH 2 OH + NH 2 CH 2 CN -> CH3OH + NHCH 2 CN 


6100 


1.0 


(#29) 


32 


CH 2 OH + NH 2 OH -> CH3OH + HNOH 


5130 


1.0 


Evans-Polanyi 


33 


CH 2 OH + N 2 H 4 -y CH3OH + NH 2 NH 


6100 


1.0 


(#29) 


34 


CH 2 OH + NH2OCH3 -> CH3OH + HNOCH3 


6100 


1.0 


(#29) 


35 


CH 2 OH + NH 2 CH 2 OH ->■ CH3OH + HNCH 2 OH 


6100 


1.0 


(#29) 


36 


CH 2 OH + NH2CONH2 ->■ CH3OH + NH2CONH 


6100 


1.0 


(#29) 


37 


CH 2 OH + NH 2 COOH CH3OH + HNCOOH 


6100 


1.0 


(#29) 


38 


CH 2 OH + CH3OCONH2 ~¥ CH3OH + CH3OCONH 


6100 


1.0 


(#29) 


39 


CH 2 OH + HOCH2CONH2 -> CH3OH + HOCH2CONH 


6100 


1.0 


(#29) 


40 


CH 2 OH + NH 2 COCHO -> CH3OH + HNCOCHO 


6100 


1.0 


(#29) 


41 


CH 2 OH + CH3CONH2 -> CH3OH + CH3CONH 


6100 


1.0 


(#29) 


42 


CH 2 OH + NH 2 CH 2 COOH ->■ CH3OH + HNCH 2 COOH 


6100 


1.0 


(#29) 


43 


CH 2 OH + C2H5COOH -> CH3OH + C2H4COOH 


6100 


1.0 


(#29) 


44 


CH 2 OH + NH 2 C 2 H 5 -> CH3OH + NH2C2H4 


6100 


1.0 


(#29) 


45 


CH 2 OH + NH 2 CH 2 CHO -> CH3OH + NH 2 CH 2 CO 


3100 


1.0 


(#29) 
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Table 3 — Continued 



# 


Reaction 


Ea 


Width 


Ref. 






(K) 


(A) 




46 


CH 2 OH + NH 2 CH 2 CHO -> CH3OH + HNCH 2 CHO 


6100 


1.0 


(#29) 


47 


CH 2 OH + CH3COOH -> CH3OH + CH 2 COOH 


6100 


1.0 


(#29) 


48 


CH 3 + CH3CHO CH 4 + CH3CO 


1720 


1.0 


Evans-Polanyi 


49 


CH 3 + CH3OH -> CH 4 + CH 2 OH 


3600 


1.0 


Tsang (1987) 


50 


CH 3 + CH3OH CH 4 + CH3O 


3480 


1.0 


Tsang(1987) 


51 


CH 3 + CO -> CH3CO 


2500 


1.0 


Garrod et al. (2008) 


52 


CH 3 + H 2 CO -y CH 4 + HCO 


4440 


1.0 


Baulch et al. (1992) 


53 


CH 3 + HCOOCH3 -> CH 4 + CH3OCO 


5040 


1.0 


Good & Francisco (2002) 


54 


CH 3 + HCOOCH3 -> CH 4 + HCOOCH2 


5800 


1.0 


Good & Francisco (2002) 


55 


CH 3 + HCOOH -> CH 4 + COOH 


4550 


1.0 


Evans-Polanyi 


56 


CH 3 + NH 2 CHO CH 4 + NH 2 CO 


3520 


1.0 


Boden & Back (1970) 


57 


CH 3 + CH 2 (OH)CHO CH 4 + CH 2 (OH)CO 


1880 


1.0 


Evans-Polanyi 


58 


CH 3 + HNCO CH 4 + OCN 


1200 


1.0 


Estimate 


59 


CH3O + CH3CHO -> CH3OH + CH3CO 


1360 


1.0 


Evans-Polanyi 


60 


CH3O + CO -> CH3OCO 


2500 


1.0 


Garrod et al. (2008) 


61 


CH3O + H 2 CO -> CH3OH + HCO 


1500 


1.0 


Tsang & Hampson (1986) 


62 


CH3O + HCOOCH3 -> CH3OH + CH3OCO 


3200 


1.0 


Evans-Polanyi 


63 


CH3O + HCOOCH3 CH3OH + HCOOCH2 


4120 


1.0 


Thynne & Gray (1963) 


64 


CH3O + HCOOH ->■ CH3OH + COOH 


3170 


1.0 


Evans-Polanyi 


65 


CH3O + NH 2 CHO ->■ CH3OH + NH 2 CO 


1980 


1.0 


Evans-Polanyi 


66 


CH3O + CH 2 (OH)CHO CH3OH + CH 2 (OH)CO 


1470 


1.0 


Evans-Polanyi 


67 


CH3O + HNCO ->■ CH3OH + OCN 


1200 


1.0 


Estimate 


68 


CH3O + C 2 H 5 CHO ->■ CH3OH + C2H5CO 


1470 


1.0 


(#66) 


69 


CH3O + C 3 H 4 CH3OH + C3H3 


1100 


1.0 


Evans-Polanyi 


70 


CH3O + C 3 H 6 -> CH3OH + C 3 H 5 


965 


1.0 


Evans-Polanyi 


71 


CH3O + CH 2 NH -> CH3OH + H 2 CN 


1200 


1.0 


(#67) 


72 


CH3O + CH3NH2 CH3OH + CH 2 NH 2 


2530 


1.0 


Evans-Polanyi 


73 


CH3O + CH3CN -> CH3OH + CH 2 CN 


2070 


1.0 


Evans-Polanyi 


74 


CH3O + C 2 H 5 CN CH3OH + C 2 H 4 CN 


2000 


1.0 


Evans-Polanyi 


75 


CH3O + C 3 H 7 CN -> CH3OH + C 3 H 6 CN 


2000 


1.0 


(#74) 


76 


CH3O + NH 2 CH 2 CN CH3OH + NHCH 2 CN 


2000 


1.0 


(#74) 


77 


CH3O + NH 2 OH -> CH3OH + HNOH 


2000 


1.0 


(#74) 


78 


CH3O + N 2 H 4 ->■ CH3OH + NH 2 NH 


2000 


1.0 


(#74) 


79 


CH3O + NH2OCH3 -> CH3OH + HNOCH3 


2000 


1.0 


(#74) 


80 


CH3O + NH 2 CH 2 OH -> CH3OH + HNCH 2 OH 


2000 


1.0 


(#74) 


81 


CH3O + NH2CONH2 ->■ CH3OH + NH2CONH 


2000 


1.0 


(#74) 


82 


CH3O + NH 2 COOH CH3OH + HNCOOH 


2000 


1.0 


(#74) 


83 


CH3O + CH3OCONH2 -> CH3OH + CH3OCONH 


2000 


1.0 


(#74) 


84 


CH3O + HOCH2CONH2 ->■ CH3OH + HOCH2CONH 


2000 


1.0 


(#74) 


85 


CH3O + NH 2 COCHO -> CH3OH + HNCOCHO 


2000 


1.0 


(#74) 


86 


CH3O + CH3CONH2 ->■ CH3OH + CH3CONH 


2000 


1.0 


(#74) 


87 


CH3O + NH 2 CH 2 COOH ->■ CH3OH + HNCH 2 COOH 


2000 


1.0 


(#74) 


88 


CH3O + C2H5COOH -> CH3OH + C 2 H 4 COOH 


2000 


1.0 


(#74) 


89 


CH3O + NH 2 C 2 H 5 -> CH3OH + NH 2 C 2 H 4 


2000 


1.0 


(#74) 


90 


CH3O + NH 2 CH 2 CHO -> CH3OH + NH 2 CH 2 CO 


1470 


1.0 


(#66) 
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Table 3 — Continued 



# 


Reaction 


Ea 


Width 


Ref. 






(K) 


(A) 




91 


CH3O + NH 2 CH 2 CHO -> CH3OH + HNCH2CHO 


2000 


1.0 


(#74) 


92 


CH3O + CH3COOH -»■ CH3OH + CH2COOH 


2000 


1.0 


(#74) 


93 


NH 2 + CH3CHO ->■ NH 3 + CH3CO 


1250 


1.0 


Hacket al. (1986) 


94 


NH 2 + CH3OH -y NH3 + CH 2 OH 


3280 


1.0 


Evans-Polanyi 


95 


NH 2 + CH 3 OH -y NH 3 + CH3O 


3290 


1.0 


Evans-Polanyi 


96 


NH 2 + CO -y NH 2 CO 


2500 


1.0 


Garrod et al. (2008) 


97 


NH 2 + H 2 CO -y NH 3 + HCO 


2360 


1.0 


Evans-Polanyi 


98 


NH 2 + HCOOCH3 -y NH 3 + CH3OCO 


3360 


1.0 


Evans-Polanyi 


99 


NH 2 + HCOOCH 3 -y NH 3 + HCOOCH 2 


2690 


1.0 


(#113) 


100 


NH 2 + HCOOH -y NH 3 + COOH 


3340 


1.0 


Evans-Polanyi 


101 


NH 2 + NH 2 CHO -y NH 3 + NH 2 CO 


2670 


1.0 


Evans-Polanyi 


102 


NH 2 + CH 2 (OH)CHO -y NH 3 + CH 2 (OH)CO 


2400 


1.0 


Evans-Polanyi 


103 


NH 2 + HNCO -y NH 3 + OCN 


1200 


1.0 


Tomeczek & Gradon (2003) 


104 


NH 2 + CH 4 -y NH 3 + CH 3 


4680 


1.0 


Moller & Wagner (1984) 


105 


NH 2 + C2H5CHO -y NH3 + C2H5CO 


2400 


1.0 


(#102) 


106 


NH 2 + C 2 H 6 -»■ NH 3 + C 2 H 5 


5340 


1.0 


Ehbrecht et al. (1987) 


107 


NH 2 + C3H4 -y NH 3 + C3H3 


2200 


1.0 


Evans-Polanyi 


108 


NH 2 + C 3 H 6 -y NH3 + C 3 H 5 


2160 


1.0 


Khe & Lesclaux (1979) 


109 


NH 2 + C 3 H8 -y NH 3 + C 3 H 7 


4730 


1.0 


Ehbrecht et al. (1987) 


110 


NH 2 + CH 2 NH -y NH3 + H 2 CN 


1200 


1.0 


(#103) 


111 


NH 2 + CH3NH2 -y NH 3 + CH 2 NH 2 


2990 


1.0 


Evans-Polanyi 


112 


NH 2 + CH 3 CN -y NH 3 + CH 2 CN 


2730 


1.0 


Evans-Polanyi 


113 


NH 2 + C2H5CN -y NH3 + C2H4CN 


2690 


1.0 


Evans-Polanyi 


114 


NH 2 + C3H7CN -y NH3 + C 3 H 6 CN 


2690 


1.0 


(#113) 


115 


NH 2 + NH 2 CH 2 CN -y NH 3 + NHCH 2 CN 


2690 


1.0 


(#113) 


116 


NH 2 + NH 2 OH -y NH3 + HNOH 


2690 


1.0 


(#113) 


117 


NH 2 + N 2 H 4 -y NH 3 + NH 2 NH 


2690 


1.0 


(#113) 


118 


NH 2 + NH2OCH3 -y NH 3 + HNOCH 3 


2690 


1.0 


(#113) 


119 


NH 2 + NH 2 CH 2 OH -»■ NH 3 + HNCH 2 OH 


2690 


1.0 


(#113) 


120 


NH 2 + NH 2 CONH 2 -y NH 3 + NH 2 CONH 


2690 


1.0 


(#113) 


121 


NH 2 + NH2COOH -y NH3 + HNCOOH 


2690 


1.0 


(#113) 


122 


NH 2 + CH 3 OCONH 2 -y NH 3 + CH 3 OCONH 


2690 


1.0 


(#113) 


123 


NH 2 + HOCH2CONH2 -y NH 3 + HOCH2CONH 


2690 


1.0 


(#113) 


124 


NH 2 + NH 2 COCHO -y NH 3 + HNCOCHO 


2690 


1.0 


(#113) 


125 


NH 2 + CH3CONH2 -y NH 3 + CH3CONH 


2690 


1.0 


(#113) 


126 


NH 2 + NH 2 CH 2 COOH -y NH3 + HNCH2COOH 


2690 


1.0 


(#113) 


127 


NH 2 + C2H5COOH -y NH3 + C2H4COOH 


2690 


1.0 


(#113) 


128 


NH 2 + NH2C2H5 -y NH 3 + NH 2 C 2 H 4 


2690 


1.0 


(#113) 


129 


NH 2 + NH 2 CH 2 CHO -y NH3 + NH 2 CH 2 CO 


2400 


1.0 


(#102) 


130 


NH 2 + NH 2 CH 2 CHO -y NH 3 + HNCH 2 CHO 


2690 


1.0 


(#113) 


131 


NH 2 + CH3COOH -y NH 3 + CH2COOH 


2690 


1.0 


(#113) 


132 


NH 2 + C 2 H 5 OH -y NH 3 + C 2 H 5 


4080 


1.0 


Evans-Polanyi 


133 


O + CO -y C0 2 


1000 


1.27 


Goumans (tunnelling) 


134 


OH + CH3OH -y CH 2 OH + H 2 


359 


1.0 


Atkinson etal. (2001) 


135 


OH + CH3OH -y CH3O + H 2 


852 


1.0 


Warnatz (1984) 
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Table 3 — Continued 



# 


Reaction 


E A 


Width 


Ref. 






(K) 


(A) 




136 


OH + CH 3 CHO -y H 2 + CH3CO 


600 


1.0 


Estimate 


137 


OH + CO -> C0 2 + H 


80 


1.0 


Ruffle & Herbst (2001) 


138 


OH + CO ^ COOH 


2500 


1.0 


Garrod et al. (2008) 


139 


OH + HCOOCH3 -y H 2 + CH3OCO 


600 


1.0 


Estimate 


140 


OH + HCOOCH3 ->■ H 2 + HCOOCH2 


1000 


1.0 


Estimate 


141 


OH + HCOOH ->■ H 2 + COOH 


1000 


1.0 


Estimate 


142 


OH + NH 2 CHO -> H 2 + NH 2 CO 


600 


1.0 


Estimate 


143 


OH + NH 2 CHO -»■ H 2 + HNCHO 


591 


1.0 


Evans-Polanyi 


144 


OH + CH 2 (OH)CHO -y H 2 + CH 2 (OH)CO 


600 


1.0 


Estimate 


145 


OH + HNCO ->■ H 2 + OCN 


1290 


1.0 


Tsang (1992) 


146 


OH + CH 4 -y H 2 + CH 3 


1780 


1.0 


DeMore et al. (1997) 


147 


OH + NH3 -y H 2 + NH 2 


925 


1.0 


Atkinson et al. (2004) 


148 


OH + C 2 H 2 -y H 2 + C 2 H 


5280 


1.0 


Liuetal. (1988) 


149 


OH + C 2 H 4 -y H 2 + C 2 H 3 


2990 


1.0 


Baulch et al. (1992) 


150 


OH + C 2 H 5 CHO -y H 2 + C 2 H 5 CO 


600 


1.0 


Estimate 


151 


OH + C 2 H 6 -y H 2 + C 2 H 5 


1000 


1.0 


Atkinson et al. (2001) 


152 


OH + C3H4 -y H 2 + C3H3 


650 


1.0 


Warnatz (1984) 


153 


OH + C 3 H 6 -y H 2 + C 3 H 5 


730 


1.0 


Tsang (1991) 


154 


OH + C3H8 -y H 2 + C3H7 


1310 


1.0 


Hu et al. (1997) 


155 


OH + CH 2 NH -y H 2 + H 2 CN 


591 


1.0 


(#143) 


156 


OH + CH 3 NH 2 -y H 2 + CH 2 NH 2 


1000 


1.0 


Estimate 


157 


OH + CH3CN -y H 2 + CH 2 CN 


1030 


1.0 


Kurylo & Knable (1984) 


158 


OH + C 2 H 3 CN -y H 2 + C 2 H 2 CN 


4000 


1.0 


Evans-Polanyi 


159 


OH + C 2 H 5 CN -y H 2 + C 2 H 4 CN 


1000 


1.0 


Estimate 


160 


OH + C3H7CN -y H 2 + C 3 H 6 CN 


1000 


1.0 


Estimate 


161 


OH + NH 2 CH 2 CN -y H 2 + NHCH 2 CN 


591 


1.0 


(#143) 


162 


OH + NH 2 OH -y H 2 + HNOH 


2140 


1.0 


Espinosa-Garcia et al. (1993) 


163 


OH + N 2 H 4 -y H 2 + NH 2 NH 


591 


1.0 


(#143) 


164 


OH + NH 2 OCH 3 -y H 2 + HNOCH3 


591 


1.0 


(#143) 


165 


OH 4- NT-TiTT-r^nH — v T4iO 4- HNPT4tD14 


591 


1 


(#143) 


166 


OH + NH 2 CONH 2 -y H 2 + NH 2 CONH 


591 


1.0 


(#143) 


167 


OH + NH 2 COOH -y H 2 + HNCOOH 


591 


1.0 


(#143) 


168 


OH + CH 3 OCONH 2 -y H 2 + CH3OCONH 


591 


1.0 


(#143) 


169 


OH + HOCH 2 CONH 2 -y H 2 + HOCH 2 CONH 


591 


1.0 


(#143) 


170 


OH + NH 2 COCHO -y H 2 + HNCOCHO 


591 


1.0 


(#143) 


171 


OH + CH3CONH2 -y H 2 + CH3CONH 


591 


1.0 


(#143) 


172 


OH + NH 2 CH 2 COOH -y H 2 + HNCH 2 COOH 


591 


1.0 


(#143) 


173 


OH + CH 3 CH 2 COOH -y H 2 + C 2 H 4 COOH 


1000 


1.0 


Estimate 


174 


OH + NH 2 C 2 H 5 -y H 2 + NH 2 C 2 H 4 


1000 


1.0 


Estimate 


175 


OH + NH 2 CH 2 CHO -y H 2 + NH 2 CH 2 CO 


600 


1.0 


(#136) 


176 


OH + NH 2 CH 2 CHO -y H 2 + HNCH 2 CHO 


591 


1.0 


(#143) 


177 


OH + CH3COOH -y H 2 + CH 2 COOH 


1000 


1.0 


Estimate 


178 


OH + C 2 H 5 OH -> H 2 + C 2 H s O 


1510 


1.0 


Evans-Polanyi 



Note. — Reactions marked Evans-Polanyi use the parameters provided by Dean & Bozzelli (1999) unless 
otherwise marked with the reference reactions used to fit the Evans-Polanyi relationship. Reactions for which 
insufficient data exist for an Evans-Polanyi are assigned the activation-energy value of the similar reaction indi- 
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cated. Reactions marked estimate are assigned a generic value appropriate to that reaction type. 
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Table 4. Binding energies and diffusion barriers used for a selection of important atoms and 
molecules, representative of adhesion to an amorphous ice surface. 



Species 






£dif 


u 
1 1 


4^0 


^1 S 


1 ^8 

1JO 


tl 2 


HOU 


JU 1 


1 SO 


V 7 1 1 


ZOJu 


1QQ^ 

177J 


QQ8 

yyo 


n 2 w 


J / uv 


J77U 


1QQ^ 

Lyy j 


v 1 1 2 


luJU 


7^5 

/ j j 


%8 


fH, 
*- n i 


1 175 

11/ J 


893 


41 1 

Til 


v--rl4 


1 ^00 

1 jUU 


Q1 


4-JJ 


.Nil 




1UUU 


UJJ 






9779 


1 386 




55^0 


^87 1 

JO / 1 


1"JU 


V I. ) 


I 1 ^0 

I I jU 


oUJ 




CO-> 
v - v ' _ 


2575 




901 

y \j 1 


HCO 


1600 


1120 


560 


I 1 2 V V 7 




14^ 


718 




9500 


1750 

1 / J \J 


875 

O / J 


V- 1 12 Wl 1 


5080 


3556 

JJJU 


1778 

1 / / O 


v 1 1 ; v -* 1 1 




^871 
JO 1 1 


1Q^6 

1 7JU 




9900 


1 540 

1 J*TV/ 


770 


V^l 1 ^ V- v.7 


9395 

Z.JZ..J 


1 698 


814 

O It 


V- 1 1 \ V- 1 1 v.7 


9775 


1943 


971 

y 1 1 


Horo 


51 90 

J 1 Z.U 


3584 
j j 0^ 


1799 


n v U \J n 


J J /U 


^8QQ 

J07 7 


1 cKn 


CH3OCH3 


3675 


2573 


1286 


C 2 H 5 OH 


6260 


4382 


2191 


HCOOCH3 


5200 


3640 


1820 


CH 2 (OH)CHO 


6680 


4676 


2338 


CH3COOH 


6300 


4410 


2205 


(CH 3 0) 2 


5000 


3500 


1750 


CH3OCH2OH 


7580 


5306 


2653 


(CH 2 OH) 2 


10200 


7140 


3570 


HCOOC 2 H 5 


6250 


4375 


2188 


C 2 H 5 COOH 


9060 


6342 


3171 


C 2 H 5 CHO 


5540 


3878 


1939 
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Table 4 — Continued 



Species 




^swap 


Edit 


CH 3 CN 


4680 


3276 


1638 


C2H3CN 


4640 


3248 


1624 


C2H5CN 


5540 


3878 


1939 


C3H7CN 


7240 


5068 


2534 


CH 2 NH 2 


4680 


3276 


1638 


CH3NH2 


5130 


3591 


1796 


NH 2 CHO 


5560 


3892 


1946 


NH 2 OH 


6810 


4767 


2384 


(NH 2 ) 2 


7910 


5537 


2769 


NH2OCH3 


6460 


4522 


2261 


NH 2 CH 2 OH 


9040 


6328 


3164 


NH 2 CH 2 CN 


5480 


3836 


1918 


NH 2 CH 2 COOH 


10100 


7070 


3535 


NH 2 C 2 H 5 


7890 


5523 


2762 


NH 2 CH 2 CHO 


6610 


4627 


2314 


(NH 2 ) 2 CO 


7910 


5537 


2769 


NH2COOH 


9080 


6356 


3178 


NH 2 COCHO 


6710 


4697 


2349 


CH3OCONH2 


7610 


5327 


2664 






71 33 


3567 


CH3CONH2 


6280 


4396 


2198 


CH3OCOOH 


7620 


5334 


2667 


CH 2 (OH)COOH 


10200 


7140 


3570 


CH3OOH 


5350 


3745 


1873 


CH 2 (OH) 2 


7930 


5551 


2776 


(CH 3 ) 2 CO 


3500 


2450 


1225 


CH3OCOCH3 


4830 


3381 


1691 


CH3COCH2OH 


7410 


5187 


2594 


(CH 3 ) 2 CO 


6150 


4305 


2153 


CH3OCOCH2OH 


8730 


6111 


3056 


(CH 2 OH) 2 CO 


11300 


7910 


3955 
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Table 5. A selection of gas-phase reactions relevant to glycine included in the chemical network 



Reaction Rate 



NH 2 CH 2 C00H + hv (CR) 


-> 


NH 2 CH 2 CO + OH 


1.50 x 


10 3 Co 


NH 2 CH 2 C00H + hv (CR) 


-> 


CH 2 NH 2 + HOCO 


5.00 x 


10 : 


Co 


NH 2 CH 2 C00H + hv (CR) 


— » 


CH 2 COOH + NH 2 


1.00 x 


HT 


Co 


NH 2 CH 2 C00H + hv (CR) 


— > 


HNCH 2 COOH + H 


1.00 x 


10 3 Co 


NH 2 CH 2 C00H + hv 


-> 


NH 2 CH 2 CO + OH 


5.00 x 


10" 


10 exp(-1.70xA v ,)s- 1 


NH 2 CH 2 C00H + hv 


-> 


CH 2 NH 2 + HOCO 


9.00 x 


10 - 


10 exp(-1.60xA v ,) S - 1 


NH 2 CH 2 C00H + hv 


-¥ 


CH 2 COOH + NH 2 


1.00 x 


10-° 9 exp(-1.70xA v )s- 1 


NH 2 CH 2 C00H + hv 


-> 


HNCH 2 COOH + H 


5.00 x 


io- 


10 exp(-1.70xA v ,)s- 1 


C + + NH 2 CH 2 COOH 


— > 


NH 2 CH 2 COOH+ + C 


1.44 x 


10" 


9 3-1 
cm s 


H 3 + + NH 2 CH 2 COOH 


— ► 


NH 2 CH 2 COOH 2 + + H 2 


2.72 X 


10" 


9 3-1 
cm s 


H 3 0+ +NH 2 CH 2 COOH 


— > 


NH 2 CH 2 COOH 2 + + H 2 


1.19 x 


10" 


9 3-1 
cm s 


HCO+ + NH 2 CH 2 COOH 


-> 


NH 2 CH 2 COOH 2 + + CO 


1.01 x 


10" 


9 cm 3 s _1 


He+ + NH 2 CH 2 COOH 


— > 


CH 2 COOH+ + NH 2 + He 


3.97 x 


10" 


10 cmh- 1 


He+ + NH 2 CH 2 COOH 




NH 2 CH 2 CO + OH+ + He 


3.97 x 


10" 


10 cmh- 1 


He+ + NH 2 CH 2 COOH 


— > 


CH 3 NH+ + HOCO + He 


3.97 x 


10" 


10 cmV 1 


He+ + NH 2 CH 2 COOH 


— » 


CH 2 COOH + NH 2 + + He 


3.97 x 


10" 


10 cmh- 1 


He+ + NH 2 CH 2 COOH 




NH 2 CH 2 CO+ + OH + He 


3.97 x 


10" 


10 cmV 1 


He+ + NH 2 CH 2 COOH 




CH 2 NH 2 + COOH+ + He 


3.97 x 


10" 


10 cm 3 s- 1 


NH 2 CH 2 COOH 2 + + e" 


— > 


NH 2 + CH 2 + C0 2 + H 2 


1.05 x 


10" 


7 (r/300 K)-° 5 cm 3 s- 


NH 2 CH 2 COOH 2 + + e" 


— > 


NH 2 + CH 2 + HOCO + H 


3.00 x 


nr 


8 (r/300 K)-° 5 cm 3 s- 


NH 2 CH 2 COOH 2 + + e - 


— > 


CH 2 NH 2 + CO + OH + H 


3.00 X 


10" 


8 (7/300 K)-°- 5 cm 3 s- 


NH 2 CH 2 COOH 2 + + e~ 


— > 


CH 2 NH 2 + C0 2 + H 2 


3.00 X 


10" 


8 (7/300 KT - 5 cm 3 s" 


NH 2 CH 2 COOH 2 + + e" 


— » 


CH 2 NH 2 + HOCO + H 


3.00 X 


10" 


8 (r/300 KT - 5 cmV 


NH 2 CH 2 COOH 2 + + e- 


— > 


CH 2 COOH + NH 2 + H 


3.00 X 


io- 


8 (r/300 K)-°- 5 cm 3 s- 


NH 2 CH 2 COOH 2 + + e- 


— > 


NH 2 CH 2 CO + H 2 


1.50 X 


10" 


8 (r/300 KT - 5 cmV 


NH 2 CH 2 COOH 2 + + e~ 




CH 2 NH 2 + HCOOH 


1.50 X 


i(r 


8 (r/300 KT - 5 cmV 


NH 2 CH 2 COOH 2 + + e - 


-> 


NH 2 CH 2 COOH + H 


1.50 X 


io- 


8 (r/300 KT - 5 cm 3 s- 


C + + HNCH 2 COOH 


— > 


HNCH 2 COOH+ + C 


1.44 X 


10- 


9 3-1 
cm s 


H 3 + +HNCH 2 COOH 


— ► 


NH 2 CH 2 COOH+ + H 2 


2.72 X 


10- 


9 3-1 
cm s 


H 3 0+ +HNCH 2 COOH 


— > 


NH 2 CH 2 COOH+ + H 2 


1.19x 


10- 


9 3-1 
cm s 


HCO + + HNCH 2 COOH 


— ► 


NH 2 CH 2 COOH+ + CO 


1.02 x 


10" 


9 cm 3 s - ' 


He+ + HNCH 2 COOH 


— ► 


CH 2 NH + COOH+ + He 


5.96 x 


10- 


10 cmV 1 


He+ + HNCH 2 COOH 


— > 


CH 2 COOH + NH+ + He 


5.96 x 


io- 


10 cmh- 1 


He+ + HNCH 2 COOH 


— > 


CH 2 NH+ + HOCO + He 


5.96 x 


10" 


10 cmh- 1 


He+ + HNCH 2 COOH 


— » 


CH 2 COOH+ + NH + He 


5.96 x 


ur 


10 cmV 1 


HNCH 2 COOH + hv (CR) 


— > 


CH 2 NH + HOCO 


5.00 x 


10 2 Co 


HNCH 2 COOH + hv (CR) 


-> 


CH 2 COOH + NH 


1.00 x 


10- 


Co 


HNCH 2 COOH + hv 


— ► 


CH 2 NH + HOCO 


9.00 x 


10" 


10 exp(-1.60xA F )s-' 


HNCH 2 COOH + hv 


->• 


CH 2 COOH + NH 


1.00 x 


10- 


9 exp(-1.70xAi/)s- 1 


H + HNCH 2 COOH 


— » 


CH 2 NH + HCOOH 


1.00 x 


10" 


11 3-1 
cm s 


CH 2 OH + HNCH 2 COOH 


— ^ 


NH 2 CH 2 COOH + H 2 CO 


1.00 x 


10' 


11 3-1 

11 crrrs 1 


CH 3 + HNCH 2 COOH 


— ^ 


NH 2 CH 2 COOH + H 2 CO 


1.00 x 


ur 


11 3 -1 
cm s 


HOCO + HNCH 2 COOH 


-> 


NH 2 CH 2 COOH + C0 2 


1.00 x 


10" 


11 3-1 

cm s 


HCO + HNCH 2 COOH 


->■ 


NH 2 CH 2 COOH + CO 


1.00 x 


10" 


11 3-1 

cm s 


NH 2 CH 2 COOH+ + e" 


->• 


NH 2 + CH 2 + C0 2 + H 


1.35 x 


10- 


7 (r/300 KT - 5 cmV 


NH 2 CH 2 COOH+ + e" 


-> 


NH 2 + CH 2 + CO + OH 


3.00 X 


io- 


8 (r/300 K)-° 5 cm 3 s- 
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Table 5 — Continued 



Reaction 








Rate 




NH 2 CH 2 COOH+ + e" 




CH 2 NH 2 + C0 2 + H 


3.00 X 10~ 8 


(r/300K)-° 5 


cm 3 s~ 


NH 2 CH 2 COOH+ + e" 


-> 


CH 2 NH 2 + CO + OH 


3.00 X 10 -8 


(r/300K)-° 5 


cm 3 s~ 


NH 2 CH 2 COOH + + e" 


-» 


NH 2 CO + CH 2 + OH 


3.00 x 10~ 8 


(7/300 K)-°- 5 


cm 3 s~ 


NH 2 CH 2 COOH+ + e" 


-»• 


CH 2 NH 2 + HOCO 


1.50 x 10 -8 


(r/300K)-°- 5 


cm 3 s~ 


NH 2 CH 2 COOH+ + e~ 




CH 2 COOH + NH 2 


1.50 x ur 8 


(r/300K)-°- 5 


cm 3 s~ 


NH 2 CH 2 COOH+ + e" 


— ► 


HNCH 2 COOH + H 


1.50 x 10~ 8 


(77300 Kr - 5 


cm 3 s~ 


NH 2 OH+ + CH3COOH 


-> 


NH 2 CH 2 COOH+ + H 2 


t 1.00 x 10 


-11 3-1 

cm s 




NH 3 OH+ + CH3COOH 


-»• 


NH 2 CH 2 COOH 2 + + H 2 


t 1.00 x 10 


""/l-OOx 10- 


9 cm 3 s 



Note. — The cosmic ray ionization rate, (q, is set to 1.3 X 10 s throughout the model. Ay denotes 
visual extinction. CR denotes cosmic ray induced UV radiation. 

t These reactions and rates are used only in the models presented in Section 4. 



Table 6. Warm-up phase model parameters 



Model 


Time to reach 


Time to reach 




200 K 


400 K 


Fast 


5x 10 4 


7.12 x 10 4 


Medium 


2x 10 5 


2.85 x 10 5 


Slow 


1 x 10 6 


1.43 x 10 6 
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Table 7. Selection of protostellar ice compositions from the literature, and post-collapse model 

values, as a percentage of H 2 



Species 


W33A a 


NGC 7538 IRS 9 b 


Sgr A* b 


Typical 
low-mass c 


Typical 
high-mass c 


Collapse-phase 
model values 


CO 


8 


16 


<12 


29 


13 


57 


co 2 


13 


22 


14 


29 


13 


18 


CH 4 


1.5 


2 


2 


5 


2 


1.9 


H 2 CO 


6 


4 


<3 


<2 


<2 


1.6 


CH3OH 


18 


5 


<4 


3 


4 


6.9 


NH 3 


15 


13 


20-30 


5 


5 


18 


NH 2 OH 












0/0.55 d 



a Gibb et al. (2000) 

b See Gibb et al. (2000) for original references 
c Obergetal. (2011) 

d Standard value / value using enhanced NH 2 OH production described in Sec. 4.1 
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Fig. 1. — Time-dependent fractional abundances of a selection of chemical species, produced by 
the fast warm-up timescale model. Solid lines indicate gas-phase species; dotted lines of the same 
color indicate ice-mantle (surface + bulk) abundances of the same species. 
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Fig. 2. — Time-dependent fractional abundances of a selection of chemical species, produced by 
the medium warm-up timescale model. Solid lines indicate gas-phase species; dotted lines of the 
same color indicate ice-mantle (surface + bulk) abundances of the same species. 
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Fig. 3. — Time-dependent fractional abundances of a selection of chemical species, produced by 
the slow warm-up timescale model. Solid lines indicate gas-phase species; dotted lines of the same 
color indicate ice-mantle (surface + bulk) abundances of the same species. 
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T (K) 

20 50 100 200 400 



-6 

I 

C 

\ 

2 -10 
cn 
o 

-12 



-14 

5x10" 1x10 5 2x10 5 2.85x10 5 

Time (yr) 

Fig. 4. — Time-dependent gas-phase and ice-mantle fractional abundances of hydroxylamine 
(NH 2 OH), acetic acid (CH 3 COOH) and glycine (NH 2 CH 2 COOH), using the medium warm-up 
timescale model with the addition of a fast gas-phase glycine-formation mechanism and a high ini- 
tial NH 2 OH abundance in the ice mantles, as discussed in Sec. 4.1. Solid lines indicate gas-phase 
species; dotted lines of the same color indicate ice-mantle (surface + bulk) abundances of the same 
species. 
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Fig. 5. — Simulated emission spectra for the hot-core source NGC 6334 IRS 1 at 241 .350 - 241 .425 
GHz. Upper panel: Local, unconvolved emission for the on-source position (i.e. assuming pencil 
beam/infinite spatial resolution). Middle panel: Convolved emission from entire source, assuming 
a beam width encompassing the strongest glycine-emission region and appropriate to the ALMA 
Telescope (0.4 arcsec). Lower panel: Convolved emission from entire source, assuming a beam 
width appropriate to the James Clerk Maxwell Telescope (~ 20.3 arcsec); note different temper- 
ature scale of this panel. All spectra represent on-source pointings. Glycine emission at 241.373 
GHz is highlighted in red. 



2000 4000 6000 8000 10000 

r (AU) 



2000 4000 6000 8000 10000 

r (AU) 



Fig. 6. — Predicted molecular line-strength profiles for lines of sight toward NGC 6334 IRS 1 at 
radial offsets r, assuming a pencil beam. Left panel: solid line shows methyl formate emission 
at 221.979 GHz (dotted line shows glycine emission of right panel, for direct comparison); Right 
panel: glycine emission at 218.105 GHz. Line intensity is expressed as brightness temperature. 




-2-1012 -2-1012 
x (orcsec) x (orcsec) 



Fig. 7. — Predicted molecular line-strength maps for lines of sight toward NGC 6334 IRS 1, 
assuming a pencil beam (i.e. unconvolved). Left panel: methyl formate emission at 221.979 
GHz; Right panel: glycine emission at 218.105 GHz. Line intensity is expressed as brightness 
temperature; temperature scales are different in each panel. 
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Table 8. Peak fractional abundances and corresponding temperatures for a selection of 

molecules, for each warm-up model 



Fast Medium Slow 



Molecule 


n[i]/n H 


T 

(K) 


n[l]/«H 


T 

(K) 


n[i]/n H 


T 

(K) 


H 2 


1.6e-4 


392 


1.6e-4 


395 


1.6e-4 


398 


CO 


9.5e-5 


400 


1.0e-4 


391 


1.2e-4 


398 


CO2 


2.9e-5 


58 


3.0e-5 


56 


3.5e-5 


55 


CH 4 


4.4e-6 


400 


4.4e-6 


398 


3.1e-6 


27 


HCO + 


7.1e-10 


36 


7.7e-10 


96 


8.7e-10 


56 


H 2 CO 


4.6e-6 


43 


3.7e-6 


42 


5.5e-7 


39 


CH3OH 


l.le-5 


123 


9.1e-6 


122 


3.8e-6 


121 


NH3 


2.7e-5 


124 


2.7e-5 


122 


2.7e-5 


122 


H 2 2 


9.0e-7 


127 


7.6e-7 


123 


3.7e-7 


122 


OH 


4.7e-8 


70 


1 .4e-7 


50 


2.9e-8 


40 


2 


2.2e-7 


10 


5.4e-7 


109 


5.5e-6 


115 


O3 


8.8e-9 


10 


8.8e-9 


10 


8.8e-9 


10 


NO 


2.2e-6 


97 


2.0e-6 


54 


1.5e-6 


48 


HNO 


6.1e-7 


42 


4.0e-7 


41 


1.7e-7 


39 


HNCO 


8.5e-9 


400 


1.2e-7 


238 


l.le-8 


222 


OCN 


2.0e-8 


108 


9.4e-8 


59 


8.9e-9 


40 


SO 


5.4e-9 


400 


2.9e-8 


395 


3.8e-8 


400 


so 2 


7.4e-8 


76 


7.6e-8 


81 


7.8e-8 


127 


C 2 H 


1.4e-9 


27 


1.9e-9 


26 


2.5e-9 


26 


C 2 H 2 


4.5e-8 


400 


6.9e-8 


53 


7.4e-8 


205 


C 2 H 4 


2.5e-8 


71 


1.8e-8 


398 


1.2e-7 


156 


C 2 H 6 


4.6e-8 


94 


1.6e-7 


93 


9.2e-7 


90 


C 3 H 6 


1.7e-9 


400 


1.4e-8 


398 


6.6e-8 


285 


C 3 H 8 


l.le-7 


132 


1.6e-7 


130 


3.4e-7 


129 


CN 


2.6e-9 


103 


1.2e-8 


46 


4.3e-8 


41 


HCN 


3.9e-7 


338 


3.3e-7 


398 


8.2e-7 


400 


HNC 


1.5e-7 


400 


1.9e-7 


398 


7.1e-7 


400 


H 2 CN 


2.1e-10 


65 


1.6e-ll 


47 


l.le-11 


27 


CH3CN 


4.9e-9 


400 


2.1e-9 


92 


1.3e-9 


90 


C 2 H 3 CN 


l.le-8 


108 


5.6e-9 


398 


5.6e-9 


202 


C 2 H 5 CN 


5.7e-9 


123 


7.7e-8 


122 


3.0e-8 


122 


C3H7CN 


6.0e-9 


153 


1.0e-8 


149 


8.4e-10 


144 


CH 2 CO 


9.0e-8 


78 


1.9e-8 


77 


l.le-8 


203 


CH3CHO 


3.1e-9 


400 


9.4e-9 


337 


1.4e-8 


154 


C 2 H 5 CHO 


6.9e-9 


123 


1.2e-8 


120 


2.6e-8 


116 


CH 2 NH 


l.le-8 


400 


2.4e-8 


398 


1.5e-8 


161 


CH 3 NH 2 


8.0e-8 


114 


3.6e-8 


105 


6.8e-9 


100 


NH 2 CHO 


3.9e-7 


123 


1.5e-7 


122 


2.9e-9 


117 


NH 2 CN 


2.3e-8 


267 


3.6e-8 


205 


3.3e-8 


195 


NH 2 OH 


1.6e-10 


145 


6.0e-10 


141 


1.7e-9 


135 


C 2 H 5 OH 


5.9e-8 


135 


5.8e-8 


131 


4.9e-8 


127 


CH3OCH3 


4.8e-8 


400 


6.3e-8 


301 


2.8e-8 


143 
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Table 8 — Continued 





Fast 




Medium 




Slow 




Molecule 


n[i]/nu 


T 


n[i\/nn 


T 


n[i]/n u 


T 






(K) 




(K) 




(K) 




y.ze-o 


1 9ft 
1ZU 


Q 7 P 

j. /e-o 


1 1 1 

111 


a ip o 
j. ie-y 


1 03 

IUj 


iiv a jy. jv 2 ^15 


j.4e-y 


1 3A 
1 JO 


A 1 p 1 ft 

4. ic-iu 


1 31 

1 J 1 


4 ip 19 

4. ie- iz 


1 99 
1ZZ 


riv. .wwri 


A 8p 8 


1 93 


A ftp 8 

4.oe-o 


1 99 
1ZZ 


A 3p 8 
4. jc-o 


1 91 
1Z1 


(. n;,v. UUn 


1 Op in 
1 .ue-iu 


1 3zl 
1 J4 


s ftp 1 n 
j.ue- iu 


1 11 

1 J 1 


7 7p n 

/ .ze-y 


1 3ft 
1 JyJ 


r,nrnni-i 
l 2 n 5 a j i, j n 


1 . le-y 




9 7p 1 ft 

z. /e-iu 


1 on 
iyu 


3. je-y 


1 83 


I 11 AAA rl W 


j.4e-o 


Si 1 
1 


3 Dp 8 


7A 

/o 


3 ^p 

j. je-y 


7A 
/ 4 




1 5e-9 


72 


1 0e-8 


68 


1 Of 
1 .ye-y 


AA 
OO 


/" ■< IT pnMU 


z.ise-o 


1 J4 


7 7q y 
Z.Ze-is 


1 1A 


/l 7a 

4.ze-o 


1 7< 

IZj 


c^lj , orrir'T-r^ 
L,rl3 UL.UL.rl 3 


1 9p 

i.ze-y 


1 ft9 
1UZ 


7 Op 1 1 

/.ye-i 1 


0/1 

y4 


1 Hp 1 n 

i.ue-iu 


09 

yz 


^ n 2 1 1 -' 11 JLUL11.3 


8p r 


1 SO 


o.oe-y 


1 sa 

1 J4 


z.je- iz 


1 Al 
14 / 


nULUUn 


4.ze-i j 


1 71 


1 1 1 1 
1 . le- 1 1 


1 A7 
10/ 


1 ip in 
1 . ie- iu 


1 A1 
101 


HPOPOPT-TO 


1 2e-8 


Q1 


5 7e-9 


8S 

J 


2 7e-9 


82 


Ulil A nU 


a io 7 


AO 


1 3« 7 
j. je- / 


A/1 


A Ip 7 

4. ->e- / 


A9 
OZ 


unLLUUn 


j. ye- 1 1 


14 / 


1 9p 7 

i.ze- / 


1 A 9 
14Z 


9 Op 

z.ze-y 


1 38 


i\ n a la riij 


A Ap 8 

o.oe-o 


IAS 
It- J 


9 ^p 
z.oe-o 


141 
141 


1 4p 8 


1 3A 


IN HA A A J r1 


i.oe-iu 




1 9p q 

i.ze-y 


1 88 
I So 


n Op in 

y.oe-iu 


1 89 
loZ 




7 ip 


1 AO 

loy 


1 9p 7 

1 .ze- / 


1 A3 
IOj 


1 9p 8 

j.ze-o 


1 S7 

1 J / 


mu a ru^ ptro 
in n.2 v n 2 l. n w 


/ . je-y 


141 


1 2e-8 


1 3A 
1 JO 


9 Qp 
z.ye-y 


134 


IN rl2 v H 2 LIN 


A 3p Q 

4.je-y 


1 93 


1 Ip- 8 

1 .ze-o 


1 99 

1ZZ 


A Ip- 

4. /e-y 


1 1 8 

1 1 


in n.2 v n 2 v a / v / n 


8 4e- 1 1 


91 A 

ZIO 


4 6e-10 


907 
zo / 


8 le-9 


203 


lNM2L.2rl5 


A ftp Q 

o.ue-y 


1 7ft 
1 /U 


1 (la 

1 .ue-y 


1 AA 
104 


a Qp in 
o.ye-iu 


1 An 

10U 


(NH 2 ) 2 


1.6e-8 


172 


3.2e-8 


165 


3.3e-8 


160 


CH3ONH2 


2.2e-8 


140 


1.5e-8 


136 


4. le-9 


131 


CH 2 (OH)NH 2 


5.5e-8 


193 


8.5e-8 


188 


7.6e-8 


182 


CH3OCOCHO 


4.5e-ll 


118 


5.1e-13 


128 


1.6e-12 


105 


CH3OCOOH 


1.5e-12 


164 


3.8e-12 


160 


1.2e-15 


151 


CH3OCONH2 


3.2e-9 


163 


3.2e-9 


158 


1.0e-9 


153 


(CH 3 0) 2 CO 


3.4e-12 


137 


3.4e-12 


135 


1.8e-13 


127 


CH3OCOCH2OH 


1.8e-9 


187 


4.3e-10 


180 


1.0e-14 


125 


CH3OOH 


7.7e-13 


123 


3.9e-13 


120 


9.2e-14 


112 


(CH 3 0) 2 


3.2e-9 


107 


3.2e-ll 


99 


1.7e-13 


96 


CH3OCH2OH 


1.5e-8 


166 


1.4e-9 


160 


l.le-13 


151 


CH 2 (OH)CHO 


2.9e-7 


143 


l.le-7 


141 


1.9e-8 


147 


CH 2 (OH)COCHO 


3.3e-8 


240 


5.0e-8 


209 


l.le-8 


195 


CH 2 (OH)COOH 


4.5e-10 


216 


6.6e-ll 


202 


4.1e-15 


134 


CH 2 (OH)CONH 2 


1.3e-7 


215 


2.3e-7 


207 


4.9e-8 


203 


(CH 2 OH) 2 CO 


4.8e-9 


222 


2.2e-10 


203 


5.8e-13 


191 


CH 2 (OH) 2 


l.()e-9 


169 


7.6e-10 


164 


5.6e-14 


154 


(CH 2 OH) 2 


3.7e-9 


206 


l.le-9 


198 


4.1e-12 


158 
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Table 9. Simulated line intensities toward NGC 6334 IRS1, for a selection of molecules and 
beamwidths, using fast warm-up timescale model. 



Molecule 


Line 


■Supper 


Peak simulated 


Beam 


Convolved 


Convolved 


Observed 




frequency 




local intensity 


width 


intensity 


integrated intensity 


integrated intensity" 




(GHz) 


(K) 


(K) 


(arcsec) 


(K) 


(Kkms- 1 ) 


(Khns 4 ) 


Methanol 


220.079 


96.61 


113 


21.9* 


1.69 


9.03 


18 




7^0 097 


Jy.OJ 


118 


91 CP 


1 69 


9 00 


1 7 




234.683 


60.92 


116 


20.7* 


1.76 


9.40 


18 




338.125 


78.08 


107 


14.4* 


4.24 


22.6 


56 




338.345 


70.55 


106 


14.4* 


4.30 


23.0 


71 




340.141 


44.67 


120 


14.3* 


3.34 


17.8 


27 




356.007 


295.3 


107 


13.9* 


4.14 


22.1 


29 


Methyl formate 


220.190 


103.1 


246 


21.9* 


0.897 


4.78 


6.4 


(HCOOCH3) 


221.979 


154.7 


248 


21.7* 


0.980 


5.23 


5.7 




226.719 


120.2 


247 


21.3* 


1.00 


5.33 


14 




338.356 


267.2 


272 


14.4* 


1.53 


8.14 


8.8 




338.396 


257.8 


268 


14.4* 


1.63 


8.70 


11 




346.660 


320.4 


286 


14.0* 


1.20 


6.39 


7.5 




346.675 


320.4 


255 


14.0* 


1.86 


9.90 


11 




348.050 


266.1 


266 


14.0* 


1.76 


9.37 


12 


Glycine 


218.105 


201.2 


1.33 


22.0* 


2.01e-4 


0.00107 




(NH 2 CH 2 C00H) 








11. 4 C 


7.46e-4 


0.00398 












1.00 rf 


0.0869 


0.464 












0.400 d 


0.345 


1.84 






223.923 


211.9 


1.36 


21.5* 


2.12e-4 


0.00113 












11. 2 C 


7.89e-4 


0.00421 












1.00 rf 


0.0878 


0.468 












0.400 rf 


0.349 


1.86 






229.740 


223.0 


1.39 


21.0* 


2.24e-4 


0.00120 












1U.V 




U.UU44J 












1.00 rf 


0.0885 


0.472 












0.400 rf 


0.35 


1.88 






235.557 


234.3 


1.42 


20.7* 


2.33e-4 


0.00124 












10.6 C 


8.82e-4 


0.00470 












1.00 rf 


0.0891 


0.475 












0.400 rf 


0.356 


1.90 






241.373 


245.9 


1.44 


20.3* 


2.42e-4 


0.00129 












10.4 C 


9.31e-4 


0.00496 












1.00 rf 


0.0893 


0.477 












0.400 rf 


0.359 


1.91 






247.190 


257.7 


1.47 


20.0* 


2.51e-4 


0.00134 












10.1 c 


9.81e-4 


0.00523 












1.00 rf 


0.0895 


0.477 












0.400 rf 


0.360 


1.92 





Note. — Methanol line data from Xu & Lovas (1997), using partition function given by the Cologne Database for Molecular Spectroscopy 
(CDMS); methyl formate data from the JPL linelist based on the dataset of Ilyushin et al. (2009), using partition function given by JPL; 
glycine line data from Ilyushin et al. (2005), using partition function given by CDMS. 
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"Observations by Bisschop et al. (2007) 
b Appropriate to JCMT (12 m) observations at this frequency 
c Appropriate to IRAM 30 m observations at this frequency 
d Approximate value, within ALMA capabilities at this frequency 
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Table 10. Simulated line intensities toward NGC 6334 IRS1, for a selection of molecules and 
beamwidths, using medium warm-up timescale model. 



Molecule 


Line 


-Clipper 


Peak simulated 


Beam 


Convolved 


Convolved 


Observed 




frequency 




local intensity 


width 


intensity 


integrated intensity 


integrated intensity" 




(GHz) 


(K) 


(K) 


(arcsec) 


(K) 


(Kkms- 1 ) 


(Kkms -1 ) 


Methanol 


220.079 


96.61 


112 


21.9* 


1.71 


9.14 


18 




230.027 


39.83 


118 


21.0* 


1.72 


9.19 


17 




234.683 


60.92 


116 


20.7* 


1.80 


9.59 


18 




338.125 


78.08 


106 


14.4* 


4.17 


22.2 


56 




338.345 


70.55 


106 


14.4* 


4.21 


22.5 


71 




340.141 


44.67 


122 


14.3* 


3.44 


18.4 


27 




356.007 


295.3 


106 


13.9* 


4.25 


22.7 


29 


Methyl formate 


220.190 


103.1 


133 


21.9* 


0.458 


2.44 


6.4 


(HCOOCH3) 


221.979 


154.7 


166 


21.7* 


0.514 


2.74 


5.7 




226.719 


120.2 


148 


21.3* 


0.520 


2.77 


14 




338.356 


267.2 


165 


14.4* 


0.727 


3.88 


8.8 




338.396 


257.8 


170 


14.4* 


0.789 


4.21 


11 




346.660 


320.4 


148 


14.0* 


0.534 


2.85 


7.5 




346.675 


320.4 


212 


14.0* 


0.926 


4.94 


11 




348.050 


266.1 


176 


14.0* 


0.854 


4.55 


12 


Glycine 


218.105 


201.2 


2.81 


22.0* 


7.67e-4 


0.00409 




(NH 2 CH 2 COOH) 








11. 4 C 


0.00285 


0.0152 












1.00 rf 


0.321 


1.71 












0.400 rf 


1.11 


5.94 






223.923 


211.9 


2.87 


21.5* 


8.08e-4 


0.00431 












11. 2 C 


0.00301 


0.0160 












1.00'' 


0.323 


1.72 












0.400 rf 


2.10 


11.2 






229.740 


223.0 


2.93 


21.0* 


8.52e-4 


0.00454 












10 9 C 
















1.00 rf 


0.325 


1.73 












0.400 rf 


1.13 


6.05 






235.557 


234.3 


2.98 


20.7* 


8.83e-4 


0.00471 












10.6 C 


0.00334 


0.0178 












1.00 rf 


0.326 


1.74 












0.400 rf 


1.14 


6.08 






241.373 


245.9 


3.02 


20.3* 


9.14e-4 


0.00487 












10.4 C 


0.00351 


0.0187 












1.00 rf 


0.326 


1.74 












0.400 rf 


1.14 


6.10 






247.190 


257.7 


3.07 


20.0* 


9.44e-4 


0.00503 












10.1 c 


0.00369 


0.0197 












1.00 rf 


0.325 


1.73 












0.400 rf 


1.14 


6.10 





Note. — Methanol line data from Xu & Lovas (1997), using partition function given by the Cologne Database for Molecular Spectroscopy 
(CDMS); methyl formate data from the JPL linelist based on the dataset of Ilyushin et al. (2009), using partition function given by JPL; 
glycine line data from Ilyushin et al. (2005), using partition function given by CDMS. 
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"Observations by Bisschop et al. (2007) 
b Appropriate to JCMT (12 m) observations at this frequency 
c Appropriate to IRAM 30 m observations at this frequency 
d Approximate value, within ALMA capabilities at this frequency 
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Table 1 1 . Simulated line intensities toward NGC 6334 IRS 1 , for a selection of molecules and 
beamwidths, using slow warm-up timescale model. 



Molecule 


Line 


£"upper 


Peak simulated 


Beam 


Convolved 


Convolved 


Observed 




frequency 




local intensity 


width 


intensity 


integrated intensity 


integrated intensity" 




(GHz) 


(K) 


(K) 


(arcsec) 


(K) 


(Kkms- 1 ) 


(Kkms- 1 ) 


Methanol 


220.079 


96.61 


117 


21.9* 


1.65 


8.82 


18 




230 027 


39.83 


127 


21.0* 


1.63 


8.67 


17 




234.683 


60.92 


126 


20.7* 


1.73 


9.25 


18 




338.125 


78.08 


107 


14.4* 


3.97 


21.2 


56 




338.345 


70.55 


107 


14.4* 


4.00 


21.3 


71 




340.141 


44.67 


116 


14.3* 


2.99 


15.9 


27 




356.007 


295.3 


108 


13.9* 


4.16 


22.2 


29 


Methyl formate 


220.190 


103.1 


0.773 


21.9* 


0.0152 


0.0813 


6.4 


(HCOOCH3) 


221.979 


154.7 


0.792 


21.7* 


0.0161 


0.0856 


5.7 




226.719 


120.2 


0.818 


21.3* 


0.0171 


0.0910 


14 




338.356 


267.2 


0.372 


14.4* 


0.0174 


0.0929 


8.8 




338.396 


257.8 


0.416 


14.4* 


0.0194 


0.104 


11 




346.660 


320.4 


0.220 


14.0* 


0.0111 


0.0593 


7.5 




346.675 


320.4 


0.410 


14.0* 


0.0207 


0.111 


11 




348.050 


266.1 


0.417 


14.0* 


0.0208 


0.111 


12 


Glycine 


218.105 


201.2 


10.2 


22.0* 


0.00804 


0.0429 




(NH 2 CH 2 COOH) 








11. 4 e 


0.0299 


0.159 












1.00 rf 


3.20 


17.0 












0.400 rf 


8.76 


46.7 






223.923 


211.9 


10.2 


21.5* 


0.00844 


0.0450 












11. 2 C 


0.0314 


0.167 












1.00'' 


3.20 


17.1 












0.400 rf 


8.80 


46.9 






229.740 


223.0 


10.2 


21.0* 


0.00886 


0.0472 












10 9 C 


U.UJZV 


n 1 76 












1.00 rf 


3.21 


17.1 












0.400 rf 


8.82 


47.0 






235.557 


234.3 


10.2 


20.7* 


0.00914 


0.0488 












10.6 C 


0.0345 


0.184 












1.00 rf 


3.20 


17.1 












0.400 rf 


8.81 


47.0 






241.373 


245.9 


10.2 


20.3* 


0.00941 


0.0502 












10.4 C 


0.0361 


0.193 












1.00 rf 


3.18 


17.0 












0.400 rf 


8.78 


46.8 






247.190 


257.7 


10.1 


20.0* 


0.00968 


0.0516 












10.1 c 


0.0378 


0.202 












1.00 rf 


3.16 


16.9 












0.400 rf 


8.73 


46.6 





Note. — Methanol line data from Xu & Lovas (1997), using partition function given by the Cologne Database for Molecular Spectroscopy 
(CDMS); methyl formate data from the JPL linelist based on the dataset of Ilyushin et al. (2009), using partition function given by JPL; 
glycine line data from Ilyushin et al. (2005), using partition function given by CDMS. 
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"Observations by Bisschop et al. (2007) 
b Appropriate to JCMT (12 m) observations at this frequency 
c Appropriate to IRAM 30 m observations at this frequency 
d Approximate value, within ALMA capabilities at this frequency 
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Table 12. Simulated line intensities for hydroxylamine (NH 2 OH) toward W3 IRS5, using 

standard medium warm-up model 



Molecule 


Line 


^upper 


Peak simulated 


Beam 


Convolved 


Convolved 




frequency 




local intensity 


width" 


intensity 


integrated intensity 




(GHz) 


(K) 


(K) 


(arcsec) 


(K) 


(Kkms" 1 ) 


Hydroxylamine 


151.102 


46.33 


0.0354 


41.5 


4.37e-5 


4.20e-4 


(NH 2 OH) 


151.118 


14.51 


0.0391 


41.5 


5.01e-5 


4.81e-4 




164.341 


116.7 


8.05e-4 


38.2 


1.08e-6 


1.03e-5 




164.627 


94.92 


8.10e-4 


38.1 


1.12e-6 


1.07e-5 




164.883 


75.60 


7.93e-4 


38.0 


1.13e-6 


1.08e-5 



Note. — * Line data and partition function obtained from the Cologne Database for Molecular Spec- 
troscopy (CDMS). 

Appropriate to NRAO 12 m observations at this frequency. 



Table 13. Simulated line intensities for hydroxylamine (NH 2 OH) toward W3 IRS5, using 
medium warm-up model with the gas-phase glycine-formation mechanism and high initial 
NH 2 OH ice-mantle abundance discussed in Sec. 4.1. 



Molecule 


Line 


^upper 


Peak simulated 


Beam 


Convolved 


Convolved 




frequency 




local intensity 


width" 


intensity 


integrated intensity 




(GHz) 


(K) 


(K) 


(arcsec) 


(K) 


(Kkms" 1 ) 


Hydroxylamine 


151.102 


46.33 


34.7 


41.5 


0.0424 


0.407 


(NH 2 OH) 
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14.51 


37.7 


41.5 
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0.475 




164.341 


116.7 


0.883 


38.2 
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0.0107 




164.627 


94.92 


0.885 


38.1 


0.00115 


0.0111 




164.883 


75.60 


0.865 


38.0 


0.00117 


0.0112 



Note. — * Line data and partition function obtained from the Cologne Database for Molecular Spec- 
troscopy (CDMS). 

Appropriate to NRAO 12 m observations at this frequency. 
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